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This  report  was  prepared  by  General  Atomic  Division,  General  Dynamics 
Corporation,  San  Diego  California,  under  Contract  AF29(601)-7035.  The  research 
was  funded  by  DASA  under  Program  Element  6.16.46.01D,  Project  5710,  Subtask 
07.002,  and  by  ARPA  Order  313,  Program  Element  6. 25. 03. 01. R. 

Inclusive  dates  of  research  were  22  July  1965  to  21  July  1966.  The  report 
was  submitted  28  April  1967  by  the  AFWL  Project  Officer,  Maj  George  Spillman 
(WLRT) .  The  contractor's  report  number  is  GA-7370. 

This  final  report  on  Nuclear  Explosion  Interaction  Studies  is  being 
published  in  four  volumes.  The  volume  titles  are  as  follows:  Volume  I,  Methods 
for  Analysis  of  Radiative  Transfer;  Volume  II,  Methods  for  Analysis  of  Thermal 
Phenomena;  Volume  III,  Miscellaneous  Code  Development;  and  Volume  IV,  Phenome¬ 
nology  Studies  (classified  SECRET/RESTRICTED  DATA). 

The  first  three  volumes  are  devoted,  respectively,  to  theoretical  studies 
and  computer  code  development  in  radiative  transfer,  thermal  phenomena,  and 
miscellaneous  efforts  related  to  various  other  aspects  of  the  work.  The  fourth 
volume,  which  is  classified,  contains  the  results  of  applications  of  these 
techniques,  and  of  those  previously  developed,  to  the  study  of  fireball  growth 
and  the  interaction  of  laser  radiation  with  materials. 
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nature  of  the  existing  computer  programs  themselves,  which  continue  in  a  state 
of  development  as  long  as  they  are  in  use,  and  in  part  to  the  time  scale  involved 
in  bringing  new  programs  to  a  state  of  capability  for  solving  real  problems. 
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ABSTRACT 


A  group  of  SPUTTER  routines  has  been  modified  or  developed  to 
execute  in  a  more  nearly  accurate  and  efficient  manner  the  task  of  "rezoning" 
the  Lagrangian  zone  structure  in  a  specified  region  of  the  system  being 
analyzed.  The  automatic  rezoning  routine  AUTORZ  has  been  expanded  to 
provide  for  maintaining  a  prescribed  spatial  definition  within  a  region.  The 
control  routines  REZONE,  CONVER,  and  ZONE  have  been  modified  to  make 
use  of  the  new  routines  ZMS  and  ZRAD,  which  accomplish  conversion 
between  space  and  mass  coordinates,  and  a  special  procedure  PROFL  is 
described  which  performs  energy  interpolation  in  certain  cases.  To 
facilitate  maintenance  of  complex  computer  codes  such  as  SPUTTER,  OIL, 
and  HECTIC,  two  FORTRAN  analysis  routines  (FAR,  Programs  One  and 
Two)  have  been  developed  and  are  nearly  complete  at  this  time.  These 
routines  are  designed  for  replacement  of  FORTRAN  symbol  strings  by 
others,  and  for  compilation  of  complete  cross  reference  listings  of  variable 
and  subroutine  references..  Editing  routines  designed  for  specific  use  with 
the  data-link  terminal  facilities  now  installed  at  th1?  General  Atomic  Special 
Nuclear  Effects  Laboratory  (SNEL)  and  at  AFWL  are  also  described. 


(Distribution  Limitation  Statement  No.  2) 
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SECTION  I 

GENERATOR-REZONE  ROUTINES  FOR  SPUTTER 

In  response  to  the  needs  of  several  applications,  extensive  modifica¬ 
tions  have  been  made  in  most  of  the  SPUTTER  routines  to  establish  or 
revise  the  zone  structure  which  partitions  the  material  being  studied  into 
convenient  slices  for  computational  purposes.  The  more  fundamental  of 
these  changes  are  discussed  below. 

1.1.  AUTORZ  (Automatic  Rezone  Specification) 

1.1.1.  Introduction 

The  purpose  of  this  routine  is  to  maintain,  during  the  course  of  a 
calculation,  an  adequately  fine-grained  zone  structure  in  those  portions  of 
the  material  being  studied  which  are  of  principal  physical  interest,  and  at 
the  same  time  avoid  wasteful  use  of  small  zones  in  other  portions.  The 
routine  is  presently  adapted  to  a  restricted  class  of  problems,  namely, 
those  in  which  an  ablation  front  is  moving  into  solid  material  on  the  left, 
with  vapor  expanding  to  the  right.  Only  a  few  large  zones  are  needed  to 
describe  the  conduction  front  in  the  solid  material,  whereas  the  zone  at  the 
ablation  front  and  the  zones  in  the  vapor  for  a  certain  distance  outward 
must  be  much  smaller  if  the  objectives  of  the  calculation  are  to  be  achieved. 
The  outermost  portions  of  the  vapor  are  of  less  interest  and  can  be  de¬ 
scribed  by  a  less  detailed  zone  structure. 

Since  SPUTTER  is  a  Lagrangian  code,  the  zone  boundaries  move  with 
the  material  and  the  problem  is  particularly  serious.  Even  in  an  Eulenan 
code,  however,  a  similar  problem  arises  unless  a  mo  ing  coordifiate  system 
can  be  used  with  respect  to  which  the  ablation  front  is  stationary.  The 
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techniques  described  in  this  sectio.  may  easily  be  generalized  to  accom¬ 
modate  problems  associated  with  a  wide  variety  of  situations  other  titan 
those  associated  with  ablation;  near  shock  and  radiation  fronts,  for  example, 
zone  definition  may  be  required  which  is  difficult  or  impossible  to  maintain 
without  repeated  redefinition  of  the  zone  structure. 

AUTORZ  is  called  at  the  beginning  of  each  cycle  if  the  flag  EKLM  (19) 
is  nonzero.  The  AUTORZ  option  control,  THE't’AK  (58),  must  also  be 
nonzero  or  there  is  an  immediate  return. 

1.1.2.  The  Solid  Rezone  Procedure 

In  the  solid  material  it  is  necessary  to  set  up  a  few  zones  with  a  very 
large  mass  available  for  ablation  and  for  containing  the  conduction  front, 
and  to  split  these,  one  at  a  time,  into  small  zones  as  the  ablation  front 
advances  into  the  solid.  It  is  assumed  that  the  large  zones  initially  gener¬ 
ated  were  equal  in  mass,  although  this  is  not  strictly  essential.  The  splitting 
ratio  is  an  input  quantity,  THETAK(53),  and  is  typically  set  at  a  number  like 
10.0  or  20.0.  It  must  be  an  integer  but  has  REAL  type.  The  index  of  the 
boundary  between  the  large  zones  and  the  small  ones  formed  by  the  splitting 
operation  is  stored  as  THETAK(54);  this  is  reduced  by  unity  each  time  a 
large  zone  is  split.  It  is  usually  best  to  generate  the  problem  with  a  "lambda 
region  '  of  small  zones  already  present,  in  which  case  THETAK(54)  should 
be  initially  set  to  the  left  boundary  index  of  this  lambda  region. 

The  solid  rezone  is  executed  when  IA,  the  ablation  front  index,  has 
been  reduced  to  the  value  THETAK(54)+THETAK(57).  THETAK(57)  is  an 
input  specification,  usually  of  the  order  of  3  to  10,  which  establishes  a 
minimum  for  the  number  of  small  zones  behind  the  ablation  front.  This  is 
advisable  for  reasons  involving  the  conduction  and  source  deposition  calcula¬ 
tions  in  some  applications. 

The  Inner  Vapor  Rezone  Procedure 

In  the  analysis  of  certain  experimental  arrangements,  such  as  plasma 
probe  measurements,  good  zone  definition  must  be  maintained  within  a 
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region  of  space  defined  in  the  Eulerian  sense,  i.e.,  in  terms  of  fixed  coordi¬ 
nates  rather  than  moving  Lagrangian  zone  boundaries.  The  region  of  interest 
is  assumed  to  extend  from  the  ablation  front  on  the  left  to  a  coordinate 
position  specified  by  the  input  in  THETAK(59).  Within  this  region,  the  de¬ 
sired  zone  width  is  specified  by  THETAK(60).  If  THETAK  (59)=0  or  if  the 
leading  edge  of  the  vapor  has  not  yet  passed  beyond  THETAK(59),  no  attempt 
to  rezone  the  inner  vapor  is  made.  The  index  K  of  the  zone  which  contains 
the  outer  boundary  THETAK(59)  of  the  region  of  interest  and  the  index  J  of 
the  outermost  zone  within  this  region  whose  width  is  less  than  THETAK(60) 
are  then  identified.  Unless  the  solid  rezone  is  being  done  at  the  same  time, 
the  inner  vapor  rezone  is  omitted  when  the  width  of  zone  K  is  less  than  1.5 
times  THETAK(60).  If  the  solid  rezone  is  to  be  done  anyway,  the  inner  vapor 
rezone  is  also  done  provided  that  there  is  at  least  one  zone  between  J  and  K. 

The  inner  vapor  r.ezone  is  specified  as  follows:  left  boundary  index 
IL=J+1,  right  boundary  index  IR=K+1,  left  zone  widch  -VAL(2)  =  THETAK(60), 
number  of  zones  VAL(5)=one  more  than  the  minimum  number  for  which  none 
of  the  new  zones  would  exceed  THETAK(60)  in  width.  This  specification 
generates  new  zones  whose  width  decreases  slightly  toward  the  right,  a 
desirable  situation  since  it  compensates  in  advance  for  later  expansion  of 
the  zones. 

1.1.4.  The  Outer  Vapor  Rezone 

In  order  to  prevent  the  total  nuj^ber  of  zones  from  increasing  too 
rapidly  due  to  repeated  solid  rezoning,  the  outer  part  of  the  vapor  is  rezoned 
periodically  to  increase  zone  size  in  this  region.  The  outer  region  extends 
from  the  outer  boundary  of  the  "inner  vapor  region"  to  the  leading  edge  of 
the  vapor,  R(I3).  The  inner  vapor  region  consists  of  the  rezone  region 
defined  in  the  preceding  paragraph,  if  such  a  rezone  is  being  done;  otherwise 
it  consists  of  a  definite  (input)  number  THETAK(55)  of  zones  counted  out 
from  the  ablation  front.  The  outer  vapor  rezone  is  normally  done  whenever 
the  sclid  or  inner  vapor  rezones  are  done;  however,  if  the  inner  j apor  is 
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being  rezoned,  the  outer  region  must  have  a  minimum  of  THETAK(55)  zones 
in  it  or  it  is  not  considered  worthwhile. 

The  outer  region  is  rezoned  as  follows:  left-most  zone  mass  VAL(2) 
=the  mass  of  the  outermost  of  the  inner  zones,  ratio  of  each  zone  mass  to 
its  neighbor  on  the  left  VAL(4)=THETAK(56).  (In  order  to  reduce  the  number 
of  zones,  THETAK(56)  must  be  greater  than  1.) 

1.1. S  Additional  Controls  and  Flags 

If  THETAK(58)  =  1 only  solid  rezones  can  be  done;  if  THETAK(58)= 

-1.,  solid  rezones  are  possible  only  on  cycles  for  which  both  inner  and  outer 
vapor  rezones  are  also  done.  For  THETAK(58)=2.,  there  are  no  controls 
beyond  those  listed  in  the  preceding  paragraphs.  If  a  solid  rezone  is  to  be 
done,  the  temperature  profile  flag  is  set  (this  is  a  min.ua  sign  on  the  right 
boundary  index  IR  for  the  solid  rezone  specification).  If  any  rezones  are 
to  be  done,  the  linkage  flag  S16  is  set  to  -1,  indicating  to  the  rezone  pro¬ 
grams  that  they  are  to  rezor.o  from  data  in  RK,  and  MAIN1  is  then  called. 

1.1.6.  The  RK  Array 

The  RK  array  contains  specifications  for  a  maximum  of  nine  rezone 
regions,  which  the  rezone  program  will  process  in  the  order  of  placement 
in  the  RK  array.  Thus,  AUTORZ  specifies  that  the  solid  rezone  is  to  be 
done  after  the  vapor  rezones  by  listing  its  specifications  after  thosn  for 
the  vapor  rezones  (if  any). 

RK(1)  contains  the  number  of  rezone  regions  specified  in  the  array. 
Each  specification  consists  of  11  words,  i.e.,  RK(2)-RK(12),  RK(13)-RK(23), 
etc.  These  words  are,  respectively: 

IL  Left  boundary  index  of  region  (REAL  type) 

IR  Right  boundary  index  of  region  (REAL  type) 

VAL(2)  Mass  (if  positive)  or  width  (if  negative)  of  left-most  zone 

VAL(3)  Mass  (if  positive)  or  width  (if  negative)  of  right-most 
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VAL(4)  Width  ratio  (if  VAL(6)=1)  or  ms'.s  ratio  (if  VAL(6)=2)  of 

each  zone  to  its  neighbor  on  the  left 
VAL(5)  Number  of  zones  to  be  generated  in  the  region 

VAL{6)  Mass-space  switch  (see  VAL(4)  ) 

VAL{7)  Material  identification  number 

VAL(8)  Multiplicative  factor  for  opacity 

VAL(9)  A  (mean  atomic  mass  number) 

VAL(10)  ZMEAN  (mean  atomic  number) 

Of  these  11  parameters,  those  required  are  the  firsttwo,  at  least  two  of  the 
next  four,  and  the  seventh  (VAL(6)).  The  others  can  either  be  calculated 
from  those  given  or  retain  their  previous  values  unchanged.  The  AUTORZ 
routine  always  leaves  the  last  four  words  unspecified. 

1.2.  REZONE  (REZONE  CONTROL  ROUTINE) 

On  entry  to  this  routine,  rezone  specifications  have  been  placed  in  the 
RK  array  either  by  AUTORZ  or  by  the  ZONERD  routine  which  reads  the 
specifications  from  card  input.  The  number  of  such  specifications  is  given 
by  RK(1),  and  the  entire  procedure  described  below  is  repeated  for  each 
specification  in  turn.  The  format  for  the  specifications  is  given  in  Section 

1 .1.6. 

The  rezone  procedure  will  modify  information  within  the  rezone 
region  defined  by  the  boundary  index  range  IL^I-IR.  Unmodified  values 
must  first  be  saved,  as  follows: 
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F or  all  I 

1  i 


Quantity 

Save  location 

Velocities,  RD(I) 

RDTAB(N) 

Densities,  l./SV(D 

RHOTAB(N) 

Specific  energies,  E(I) 

THETAB(N) 

Coordinates,  R(I) 

RTABRD(N) 

Coordinates,  R(I) 

RTABRO(N) 

Coordinates,  R(I) 

RTABTH(N) 

* 

Mean  ionic  charges,  FEW(I) 

CAR(N) 

Z1,  X6{I) 

CAPAC(N 

Z2,  BC(I) 

CHIC(N) 

Z3,  BR(I) 

Elemental  mean 

CHIR(N) 

Z4,  CRTC(I) 

ionic  charges 

X4(N; 

. 

Z5,  RHO(I)  . 

X5(N) 

Zone  masses, 

G(I) 

EC(N) 

Mass  depletion  energies, RFF(I) 

PB(N) 

Mass  depletion  ar«a, 

SMLQ(N) 

SMLR(I)**2 

Left  boundary  of  rezone  region, 

RR 

R(IL) 

Right  boundary  of  rezone  region, 

SAVRV6 

R(IR) 

Pusher  •  elocity  RD(1)  or 

RLF 

RD(IG-l) 

Solid-vapor  interface  (see  below) 

RPUSH 

LTE-non-LTE  interface  R(NR) 

RNR 

Left  boundary  of  active  region 

RICA 

R(ICA) 

Right  boundary  of  active  region 

RIC3 

R(ICB) 


in  /NAMLC/ 
COMMON 


in  blank 
COMMON 


/NAMEC/  or 
internal 
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The  solid-vapor  interface  is  defined  by  R(IA)  for  a  left  pusher  (solid  region) 
or  R(IB)  for  a  right  pusher,  provided  that  the  vaporization  flag  BOILA  is 
zero.  If  EOILA  t  0,  RPUSH  is  located  within  the  vaporizing  zone  at  a 
position  defined  by  the  fraction  of  the  mass  in  the  zone  which  is  still  un¬ 
vaporized.  When  the  new  zone  structure  is  established,  these  mass  frac¬ 
tions  are  recomputed  if  necessary. 

The  subroutine  GENRAT  is  then  called  to  fill  in  any  previously  un¬ 
specified  quantities  VAL(2),  (3),  (4),  (5),  and  subroutine  MOVE  is  called  to 
shift  blank  COMMON  array  elements  with  indices  I  -  IR  to  their  new 
locations.  The  re  zone  index  IR  for  the  current  rezone  region  and  the 
indices  IL  and  IR  in  the  RI<  array  for  yet  unprocessed  regions  to  the  right 
of  the  current  region  a*e  also  shifted.  Other  indices  and  quantities  which 
have  been  saved  as  described  above  are  shifted  to  their  new  values.  The 
ZONE  routine  (see  below)  generates  new  boundary  positions,  masses,  and 
specific  volumes  in  the  rezone  region.  Boundary  velocities  are  recomputed 
by  linear  spatial  interpolation. 

The  next  step  is  to  recompute,  for  the  new  zones,  the  quantities  E(I), 
FEW(I),  RDD(IV  SMLR(I),  and  the  five  elemental  mean  ionic  charges.  This 
is  done  by  mass -weighting  (except  for  RDD)  to  produce  corresponding  ex¬ 
tensive  quantities,  computing  the  amount  of  each  of  these  quantities  in  each 
subzone  (subzones  are  defined  by  the  union  of  the  old  and  the  new  zone 
boundaries),  assembling  the  subzone  quantities  belonging  to  each  of  the  new 
zones,  and  dividing  by  the  new  zone  masses  (except  for  RDD)  to  restore  the 
original  physical  dimensions.  Subroutine  UP37A  is  called  to  update  lambda- 
region  arrays  and  clear  the  "save"  locations  of  the  old  information  which  is 
no  longer  needed. 

Next,  subroutine  "PROFL"  (see  below)  is  called  if  its  flag  was  set,  to 
eliminate  by  a  special  procedure  any  plateaus  in  the  specific  energy  dis¬ 
tribution  in  the  solid  which  result  from  zone  splitting.  Finally,  for  LTE 
zones,  the  temperature  and  all  of  the  thermodynamic  variables  except  the 
specific  internal  energy  are  evaluated  by  an  inversion  of  the  state  equation. 
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If  the  material  identification  number  is  in  the  range  201  to  300,  this  requires 
an  iterative  procedure,  using  interval  halving  or  similar  techniques.  For 
non-LTE  zones  and  I.TE  zones  which  specify  the  general  ionic  state  equation 
routine  EIONX  (by  material  identification  numbers  -  200  or  >  300),  a  more 
efficient  procedure  is  being  developed  which  makes  use  of  the  established 
elemental  mean  ionic  charges  and  requires  no  iteration.  At  present,  a 
simple  version  of  this  new  procedure  is  used  for  non-LTE  zones  in  one- 
material  problems. 

1.3.  CONVER,  ZONE,  ZRAD,  AND  ZMS 

The  most  important  mathematical  problem  encountered  in  generating 
or  rezoning  a  problem  is  the  conversion  between  spatial  coordinates  and 
mass  coordinates.  The  relationship  between  these  is  specified  by  a 
density  profile,  which  is  just  a  list  of  material  densities  for  each  of  the 
original  zones  (in  the  rezoning  case),  or  is  given  by  input  in  generating  anew 
problem.  The  four  routines  described  in  this  section  are  the  ones  involving 
this  conversion  process,  and  are  therefore  discussed  together. 

Zoning  may  be  carried  out  either  with  space  or  mass  as  the  basic 
independent  parameter.  If  VAL(6),  the  mass-space  switch,  is  1.,  the  zones 
are  to  be  set  up  with  a  uniform  ratio  VAL(4)  of  spatial  widths,  while 
VAL(6)=2.  specifies  a  uniform  ratio  VAL(4)  of  zone  masses.  With  space 
zoning,  the  conversion  problem  is  one  of  integrating  the  density  profile  to 
obtain  new  zone  masses,  whereas  with  mass  zoning,  the  problem  is  to 
integrate  the  reciprocal  of  the  density  profile  to  obtain  new  zone  widths. 

The  problem  is  nontrivial  for  three  reasons: 

1.  The  space  variable  may  be  the  radius  variable  of  a  cylindrical 
or  spherical  coordinate  system,  in  which  case  the  interpretation 
to  be  given  the  density  profile  is  not  without  ambiguity. 

2.  Thi  zone  widths  must,  in  many  cases,  add  up  exactly  to  a  value 
(the  total  region  width)  prescribed  in  advance,  and  the  new 
boundary  positions  must  conform  to  the  zone  widths  as  well  as 
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possible  within  the  limits  of  the  finite  computer  word  size.  The 
sequences  of  boundary  positions  and  zone  widths  need  not  be 
strictly  derivable  from  each  other,  but  they  must  both  be  regular 
and  monotonic,  with  initial  and  terminal  terms  fixed.  In  some 
applications  the  boundary  coordinates  and  zone  widths  are  8  to 
10  orders  of  magnitude  apart,  and  numerical  methods  designed 
for  easier  problems  do  not  work. 

3.  For  reasons  of  efficiency,  neither  involved  iterative  procedures 
nor  indiscriminate  use  of  double -precision  operations  is  an 
acceptable  approach  to  these  problems. 

The  interpretation  given  the  density  profile  in  generating  a  new 
problem  is  that  the  complete  density  function  is  to  be  obtained  by  inter¬ 
polating  between  the  given  points,  and  that  the  interpolant  is  linear  not  in 
the  coordinate  r,  but  in  r  where  01  =  1,2,3  for  plane,  cylindrical,  and 
spherical  symmetry  respectively.  This  eliminates  the  need  for  solving 
anything  more  complicated  than  a  quadratic  equation  in  finding  the  width  of 
a  zone  given  the  mass.  In  rezoning,  the  density  profile  is  in  origin  a  dis¬ 
crete  mass  profile,  and  any  interpolation  at  all  would  be  in  a  sense  gratu¬ 
itous.  The  assumption  which  is  made  throughout  SPUTTER  is  that  the 
density  is  constant  throughout  each  zone;  in  rezoning,  therefore,  the  pro¬ 
file  is  considered  as  defining  a  histogram,  and  the  integration  process  is 
thus  always  linear.  The  integration  subroutines  ZMS  and  ZRAD  are  di¬ 
rected  to  these  two  alternative  interpretations  of  the  density  profile  by  a 
switch  IZSWT  which  is  set  in  MAIN1  to  -1  for  rezoning  and  to  0  for 
generating. 

The  union  of  the  set  of  new  zone  boundaries  and  the  set  of  old  zone 
boundaries  or  profile  points  thus  defines  a  set  of  subzones  within  each  of 

0t 

which  the  density  is  either  constant  or  linear  in  r  .  The  mass  integration 
routine  ZMS  is  called  for  each  subzone  to  determine  the  mass  contained 
therein.  The  inverse  problem  is  solved  by  the  ZRAD  subroutine,  in  which 
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each  subzone  has  a  given  mass  and  it  is  required  to  find  the  volume  from 
the  profile.  In  these  routines  the  evaluation  of  the  slope  from  the  profile 
is  done  in  single  precision,  while  the  integration  calculation  itself  may  be 
done  in  double  precision. 

The  CONVER  routine  is  called  from  GENRAT  to  convert  specifications 
into  consistent  form.  Thus,  if  VAL(6)=2.  specifying  mass  zoning,  but  any  of 
VAL(l),  VAL{2),  or  VAL(3)  are  specified  with  negative  sign  to  indicate  that 
they  are  widths,  not  masses,  then  the  CONVER  routine,  with  the  aid  of  ZMS 
and  the  density  profile,  computes  the  proper  masses  where  necessary  for 
VAL(l)  (region  mass),  VAL(2)  (first  zone  mass),  and  VAL(3)  (last  zone  mass). 
Similarly,  CONVER  calls  ZRAD  in  the  less  usual  case  where  space  zoning 
is  specified  in  terms  of  the  mass  of  the  region  or  of  the  first  or  last  zones. 

The  ZONE  routine  computes  the  mass,  width,  boundary  coordinates, 
and  specific  volume  of  each  zone.  After  generating  the  zones  using  VAL(4) 
as  the  generating  ratio,  the  subzone  partition  is  established  and  ZMS  (for 
space  zoning)  or  ZRAD  (for  mass  zoning)  is  called  for  each  in  turn  to  com¬ 
plete  the  description  of  the  subzone.  The  assemblage  of  the  subzones  into 
the  new  zones  to  which  they  belong  may  be  done  with  double  numerical  pre¬ 
cision.  The  largest  discrepancies  arise  for  the  boundary  positions,  which 
are  cumulative.  If  the  right  region  boundary  is  prescribed  (always  the  case 
in  rezoning  and  often  the  case  in  generating),  the  generated  region  will  ex¬ 
ceed  or  fall  short  of  the  target  by  an  amount  which  is  not  necessarily  small 
compared  with  the  width  of  the  last  zone,  even  if  double  precision  operations 
as  described  above  are  used.  Smaller  discrepancies  arise  for  the  width  or 
mass  of  the  last  zone,  since  the  prescribed  value  VAL(3)  will  not,  in  general, 
agree  to  the  last  digit  with  the  value  resulting  from  starting  with  VAL(2)  and 
multiplying  N  times  by  VAL(4). 

In  each  of  these  cases  the  discrepancies  are  removed  by  the  following 
Procrustean  device:  The  quantity  in  question  is  regenerated  in  reverse 
order  (from  right  to  left),  using  the  prescribed  right-end  value  as  the  initial 
condition.  In  calculating  widths  or  masses,  the  initial  condition  is  thus 


10 


AFWL-  TR-66  - 108,  Vol  III 

VAL(3)  and  the  generating  ratio  is  l./VAL(4),  while  in  calculating  boundary 
positions,  the  initial  value  is  R(IR)  and  the  increment  is  -DELTAR(l),  where 
DELTAR(I)  is  the  same  as  the  value  used  in  the  forward  sweep.  Each 
quantity  obtained  in  the  forward  sweep  is  then  multiplied  by  a  factor  W(I) 
and  each  reverse- sweep  quantity  by  l.-W(I),  and  the  two  results  are  then 
added  to  give  the  final  result.  W(I)  is  defined  as  that  linear  function  of  the 
index  I  which  assumes  the  value  1.  at  the  left  end  of  the  region  and  0.  at  the 
right  end. 

It  may  be  the  case  that  this  last  procedure  eliminates  the  need  for 
double  precision  arithmetic  at  any  stage  of  the  calculation.  Experimentation 
on  this  point  is  planned. 

1.4.  PROFL 

The  PROFL  subroutine  is  designed  to  eliminate  a  severe  perturbation 
of  the  solid  heat  conduction  calculation  which  results  from  splitting  a  single 
large  zone  into  10  or  20  small  zones,  each  with  the  same  specific  internal 
energy  and  therefore  the  same  temperature.  The  effect  of  this  is  to 
multiply  the  conductive  flux  at  the  ends  of  this  region  by  a  large  factor  and 
to  destroy  it  entirely  in  the  interior. 

It  is  assumed  that  the  conduction  front  is  moving  to  the  left,  so  that 
the  temperature  is  a  nondecreasing  function  of  position  r.  The  profile  to 
be  established  is  required  to  assume  the  same  values  at  the  ends  of  the 
region  as  the  adjacent  zones  outside,  and  to  conserve  the  original  total 
energy  content.  These  three  conditions  are  sufficient  to  determine  a 
quadratic  function  but  not  to  ensure  that  the  quadratic  is  a  non-decreasing 
function  of  r.  The  profile  is  therefore  defined  as  a  quadratic  which,  start¬ 
ing  £.t  the  right  boundary,  decreases  monotonically,  reaching  the  left 
boundary  value  at  a  point  which  may  be  in  the  interior.  The  profile  in  that 
case  is  continued  as  a  constant  to  the  left  boundary.  The  energy  in  each 
zone  is  then  simply  the  integral  of  this  profile  across  the  zone. 

The  geometry  may  have  plane,  cylindrical,  or  spherical  symmetry. 
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L  5-  APPENDIX,  code  listings 

subroutine  automz 

CUMPILLU  MAY  16.  1»66 

*^***^*****************M.  I,***.**,  4, 


s  P  U  T  1  L  K  (  it  m  r  t>  II 


COMMON  LMDA ( 37) #  NR 


1 

KMAX 

ULANK1 »  DLANK2 

2 

II 

IG  »  NHAD 

3 

IIP1 

IGM1  »  1  ALPHA 

4 

FREQ 

CNTMAX.  AR 

5 

CVA 

CVB  »  SLUG 

b 

EM1NA 

EM I NO  ,  CA 

7 

RHOL 

RHOR  »  LPIO 

« 

RPIA 

RPIB  »  RHUIA 

COMMON 

TD  .  TL 

1 

DTMAX1 

DTMAX2.  L.TMAX3 

2 

GAMA 

WCRIT  .  SIGMAO 

3 

ROI A 

ROI AMI »  R01D 

4 

54 

S5  ,  Sb 

5 

S12 

S13  »  S14 

6 

S20 

EO  ,  FO 

C 

c* 


ASQ 
B  DELR 
COMMON 
W 
V 

X3 

SMLA 
EC 

UIGB 
CHIR 
CRTPC 
COMMON 

1  FRLM 

2  ZP1 

3  RL 

4  HEAD 


1 

2 

3 

4 

5 

6 
7 
0 


(152). 
(  37). 

(152). 
(152)  . 
(152). 
(152). 
(152). 
(152). 
(152) , 
(152). 


37)  . 
37). 
37). 
12). 


ML) 

P 

P2 

E 

G 

X4 

SMLU 

EH 

CV 

CAPAC 

GOFR 

TELM 

WLM 

ZP2 

RHOK 

MAXL 


(152 

(152) 

(152) 

(152) 

(152) 

(152) 

(152) 

(152) 

(152) 

(152) 

(152) 

(  37) 

(  37) 

(  37) 
(104) 


IJGfll  l<  , 
i  it.  AH|\  3  t 
HLANkH . 
PLANKS, 
A5MI  It  , 
ALPHA  . 
CL)  • 
LPSI  , 
MPDIH  . 
DTI  12  . 

L)TH  » 
AC  » 
K01HP1 . 
57  , 

515  . 

TAU  . 
VD 
PI 
sv 
LI 
t> 

X5 

SMLC 
SMLO 
HC 

CAPAR 
FEW 
CKLM 
OLM 
SOLID 
RDK 
MAXLM 


I A 

I  A 1 1 1 
I  AMI 
Til 

I’ll  SNA 
I IV  A 
OA 
F’l  A 

TI'l’llJl 
DTII2P 
SWITCH 
AC  03  I 'I 
GMS 

>  S « 

>  516 

i  ZERO 
(152) , 
(152) . 
(152). 
(1521 . 
(152). 
(152) . 
(152) . 
(152). 
(152)  . 
(152)  . 
(152)  . 

(  37). 

(  37). 

(  37). 
(104)  . 


Hi 

HiP) 

1 1  |M  ) 

T  My 
I  HIM  if, 
1 1 VI  • 

611 

Kill 

TA 

Dill) 

CO 

C 1 4  VI  ?  T 
SI 

SO 
si  7 
It 


IC.A 

K  API 

HAM) 

IILAIIM, 

HOII  A 

MCA 

GL 

MO  I A 
Til 

OTPMIli 
CM  III 
SUMP  A 
S2 
SKI 
S)o 


( lb?) . 
(15?) . 
( 152) , 
U  52 )  ► 
(152) . 
(152). 
(152). 
(152). 
(152). 
(152) . 
(152), 
(152), 

(  37), 
AMASNO(  37), 
ECHCK  (  37), 
TIItTAK  (104)  , 


aid  mm  i  ii 
Aiitmiii  oi 

AllTOIKPlIi 
At  IT  Ofi  II  >0 

•  Atnouot.il 

*  AUT 000711 
♦«  AUTOOIMO 
*  ♦AUTOOO'n) 
*• AUIO0100 
t(ll  .AUTOOUO 
ICIIP)  ,  AUT 00  1 20 
K.ltMt  .AIM  00  110 
IILLIMtT.  AUTOOKtO 
HO  I L  M  .Al/TOniSO 
lien  ,  AUTO  0  It.  0 
Git  » AIM  00170 
Ml)  1 1  j  AtJTOO  lUO 

1C  AUTOOIMO 
HTI1AX  .AUT00200 
DELTA  »Atno0210 
51MI-":  ,  AUT00220 
S3  » AUT00230 
511  , AUT00240 

SI"  » AUT00250 


RDD 
Pit 
In  iO 
EK 
C 

X6 

SMLU 

SMLH 

HR 

CRTC 
CAR  . 
ELM 


C'»*******************t*****4#4##*ttt4t^^^^^^^^^^^^  ** AUTO0440 


Dl  l  1  Alt  (152)  »  AUT00260 
SMLM  ( 152) » AUTO0270 
Pill  (1521  AUT00280 
THETA  ( 152) .AUT00290 
A  (152  » AUTOO300 
X2  ( 152  I , AUT00310 
X7  ( 152) » AUT00320 
SMLF  (152) .AUT00330 
RIGA  (152) .AUT00340 
CHIC  (152) .AUT00350 
CRTH  ( 152 ) ,  AIJT00360 
OKLM  (  ?7)  AUT00370 
FCLM  (  37) , AUT00380 
CIIRNO  (  37) ,  AUTO0390 
RK  (104) .AUT00400 
TEMP  (  16) , AUT00410 
AIJT00420 
AUT00430 
** AUTO0440 
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************  INPUT  DATA  *************** 

THETAK  ( S3  USOL  ID  rezone  RATIO 
THETAK (54) =SOLID  rezone  interface 
THETAK ( 55 ) ^NUMBER  OF  SMALL  VAPOR  ZONES 
THETAK(56)=SSIGMA  IN  VAPOR  TO -RIGHT  OF  INTERFACE 
THETAK ( 57 ! SREZONE  CRITERION 
THETAK (58)  s  1.  FOR  SOLID  REZONL  ONLY, 

THETAK (58)  .EO.-l.  FOR  REZONE  IN  VAPOR 
THETAK (59)  =  RT.  BOUNDARY  FOR  CONSTANT 
THETAK (60)  =  CONSTANT  ZONE  SIZE 
IF  ( THETAK (58) .ECKO* )  GO  TO  1000 
ISRSTHETAK(54) 

NHEZsTHETAK (57) 

NSVZ  =  THETAK (55) 

IF  (IAiLE. ISR+NREZ)  GO  TO  10 
IF  (THETAK(b8) .EQ. 1. )  GO  TO  1000 
IF  (THETAK (59) .EG.O. )  GO 
IF  (R(Ib) «LT .THETAK (59) ) 

RZCF  =  1, 

GO  TO  11 
RZCF  s  0. 

THETAK ( 54 )=I SR- 1 


AUT00460 
AUTO0470 
AUTO0480 
AUT00490 
AUTO0500 
AUT00510 

„  AUT  00520 
-2.  FOR  SOLID  AND  VAPOP  AUTO0530 
ONLY  , 

ZONES 


AUTO0540 

AUTOO550 


TO 

GO 


1000 
TO  1000 


AUT00580 
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1 

I 

1 

\ 

j 


1 


II  (THLTAK(bU).Fu.l.)  GO  1*i  lim 
II-  (THLTAK(b9).r.U.O.)  GO  TO  .’ll 

11  K  =  0 

UO  12  I  :  I  At  111 
J  =  IA  ♦  IH  -  I 

II-  (R(J).GT.THLTAK(b9))  GO  10  1J 
IF  (K.EO.O)  K  =  J 

IF  (DELTAR(J).LT.TltLTAK(GU)>  GO  TO  1 A 

12  CONTINUE 

K  =  RT.  INDEX  FOR  CONSTANI  ZONI 
J  =  LF.  INDEX  FOR  CONSTANT  ZONLG 
-  -  FULL  REZONL  -  “ 

U  IF  (KZCF.GT.O.  .AND.  DLLTAR  (K )  .LT  .  1  .!<*T!!I  TAK  <  GO ) )  GO  TO  limit 
IF  (RZCF.EQ.O.  .AND.  K-J.LT. 3)  GO  TO  .’II 
RK (2)  =  J  ♦  1 
RK ( 3)  =  K  ♦  1 
RK(4)  =  -THETAK (GO) 

UK (5)  =  0. 

RK (6)  =  0. 

RK  (7)  =  INT((R(K+1)-R(J)  )  /  UH.  TAK  <011  >  >  ♦  1 
RK(fl)  =  1. 

IF  (IB.GE.IK  ♦  1  ♦  NSVZ>>  GO  TO  14 
IF  (RZCF.LO. 1 . )  GO  TO  409 
GO  TO  21 

14  RK ( 13)  =  K  ♦  1 
RK ( 14)  =  lb 
RK ( lb)  =  G (K ) 

RK ( 16)  =  0. 

RK ( 17)  =  THETAKlbG) 

RK ( lb)  =  0. 

RK ( 19)  =  2.0 

IF  IRZCF.EO.l.)  GO  TO  40 
RK (24)  =  THETAK(b4) 

RK (2b)  =  -ISR 

RK (26)  =  G(ISR-l)  /  THLTAK(b3) 

RK (27)  =  KK(26) 

RK ( 26)  =  1.0 
RK(29)  =  THETAK(b3) 

RK ( 30 )  =  2.0 
RK ( 1 )  =  3. 

GO  TO  bOO 
C  -  -  OUTER  VAPOR  FIRST  -  - 

20  IL  s  IA  +  NSVZ 

IF  (CHETAK(59).Nt.Q.)  IL  =  K 
IF  (IB.LE.IIL+NSVZ)  )  GO  TO  100 
RK (2)  =  IL 
RK(3)  =  IB 
HK(4)  =  G(IL-l) 

RK(5)  =  0. 

RK (6)  =  THETAKlbG) 

RK (7)  =  0. 

RK(fl)  =  2.0 

21  IF  (THETAK < SB) »EQ. (~1 . ) )  GO  TO  499 

C  -  SOLID  AFTER  VALOR  -  - 

30  RK (13)  =  THETAK (b4) 

HK (14)  —  —ISR 

RK(ib)  =  G( ISR-1 ) /TMtTAK I b3) 

RK ( 16)  =  RK(lb) 

RK ( 17)  =  1.0 

HK ( 18)  =  THETAK (S3) 

RK ( 19)  s  2.0 
40  HK ( 1 )  =  2. 

GO  TO  500 

C  -  -  SOLID  ONLY  -  - 

100  IF  (THETAK (56) .EO. (-1 . ) )  GO  TO  1HOO 

IF  ( IB+ INTI THETAK (53) ) .GT.  ISO)  GO  TO  200 
RK (2)  s  THETAK (54) 

HK (3)  =-ISH 


AIJT00700 

AHT00720 
AUTO0730 
AUT00740 
AUT00750 
AUT 00760 


AUT00620 


AUT0065G 
AIIT00670 
AU1 110770 


AUT  011790 
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RK<4»  =  G<  ISR“1 ) /TMETAK ( U3) 


NM5)  =  RK  !4 ) 

HK(6)  =  1. 

HK(7)  =  THETAMSJ) 

AUTOOfllO 

499 

KK(8)  =  2.0 

AUTO0B40 

KK(1)  =  1. 

500 

516  =  -1. 

CAUL  MP1 

AUTOOObO 

200 

51  =  2.0200 

AUTO0B70 

1000 

CAUL  UNCLE 

AIITOOBBO 

RETURN 

.  AUT00690 

END 

AUT00900 

f 
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c 
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c 

c 

c 

c 

c 
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SUBROUTINE  RLZONL 
COMPILED  MAY  ?6.  1966 


SUBROUTINE  REZONE  CALCULATES  NEW  ZONAL  PROPERTIES 
BETWEEN  ANY  TWO  POINTS  (IL.IR).  SHIFTS  THE  PRESENT 
ARRAYS  TO  FIT  THESE  POINTS. 

CALLED  FROM  EITHER  1MAIN  OR  AUTOZ 


RLZOOOIO 

REZ00020 

REZ00050 

REZ00060 

REZ00070 

REZOOOfiO 

REZ00090 

♦REZOOIOO 


c* 

S  P 

U  T  T  E 

R  C 

0  M  M  0 

N 

**REZ00110 

c * 

♦♦REZO0120 

c * 

**REZ00130 

COMMON 

LMDA(37>* 

NR  . 

NSMLR 

.  IA 

.  IB 

.  ICA  . 

ICB  .REZ00140 

1 

KMAX 

*  BLANK 1* 

BLANK2. 

BLANK3 

.  I API 

.  IBP1 

.  ICAP1  . 

ICBP1  .REZ00150 

2 

II 

*  IG  * 

NR  AD  . 

BLANK4 

;  I  AMI 

.  IBM1 

.  KAMI  . 

ICBM1  .REZ00160 

3 

IIP1 

*  IGM1  * 

1 ALPHA. 

BLANK5 

.  TH 

.  TMAX 

.  BLANK6. 

DELPRT .REZ00170 

4 

FREQ 

*  CNTMAX* 

AR  . 

ASMLR 

.  PUSHA 

.  PUSHB  .  BOILA  . 

BOILB  .REZQ01A0 

5 

CVA 

•  CVB  • 

SLUG  . 

ALPHA 

•  HVA 

.  HVB 

.  HCA  . 

HCB  .REZO0190 

6 

EMINA 

*  EMINB  * 

CA  . 

CB 

.  GA 

.  GB 

.  GL  . 

GR  . REZ00200 

7 

RHOL 

*  RHOR  • 

LPIO  » 

EPSI 

*  RIA 

.  RIB 

.  RDIA  . 

RDIB  .REZ00210 

a 

RPIA 

*  RPIB  * 

RPDIA  . 

RPDIB 

.  TPRINT 

.  TA 

.  TB  » 

TC  REZ00220 

COMMON 

TD  * 

TE  . 

DTH2 

.  DTH2P 

.  DTH1 

.  DTRMIN. 

DTMAX  .REZ00230 

1 

DTMAX1 *  DTMAX2  * 

DTMAX3* 

DTR 

.  SWITCH 

.  CO 

.  CMIN  . 

DELTA  .REZ00240 

2 

GAMA 

•  WCRIT  • 

SIGMAQ. 

AC 

.  AC03T4 

.  CNVRT  .  SUMRA  . 

SUMRB  .REZ00250 

3 

ROI A 

*  ROI AMI  * 

ROIB  . 

ROIBP1 

.  GMS 

.  SI 

*  S2  . 

S3  .REZ00260 

4 

S4 

*  S5  • 

56  . 

S7 

.  58 

.  S9 

.  S10  . 

Sll  .REZ00270 

5 

S12 

*  S13  • 

S14  • 

S15 

.  S16 

.  S17 

.  S18  . 

S19  . REZ00280 

6 

S20 

•  EO  • 

FO  . 

TAU 

.  ZERO 

.  R 

(152).  DELTAR( 152) .REZ00290 

7 

ASQ 

( 152) *  RD 

(152). 

VD 

( 152* . 

RDD 

(152).  SMLR 

( 152) .REZ00300 

8 

DELR 

I  37).  P 

(152). 

PI 

(152). 

PB 

(152).  PB1 

(152)  REZ00310 

COMMON 

P2 

(152). 

SV 

(152). 

RHO 

(152).  THETA 

(152) .REZO0320 

1 

W 

(152).  E 

(152). 

El 

(152). 

EK 

(152).  A 

( 152) .REZ00330 

2 

V 

(152).  G 

(152). 

D 

(152). 

C 

(152).  X2 

( 152) .KEZ00340 

3 

X3 

(152).  X4 

(152). 

X5 

(152). 

X6 

(152).  X7 

(152) .REZ00350 

4 

SMLA 

(152).  SMLB 

(152). 

SMLC 

(152). 

SMLD 

(152).  SMLE 

(152) . REZ00360 

5 

EC 

(152).  ER 

(152). 

smlq 

(152). 

SMLH 

(152).  BIGA 

(152) .REZ00370 

< 

BIGB 

(152).  CV 

(152). 

BC 

(152). 

BR 

(152).  CHIC 

( 152) .REZ00380 

7 

CHIR 

(152).  CAPAC  (152). 

CAPAR 

(152). 

CRTC 

(152).  CRTR 

(152) .REZ00390 

8 

CRTPC 

(152).  GOFR 

(152). 

FEW 

(152). 

CAR 

(152).  OKLM 

(  37)  REZ00400 

COMMON 

TELM 

( 

EKLM 

(  37). 

ELM 

(  37).  FCLM 

(  37) .REZ00410 

1 

FRLM 

(  37).  WLM 

(  37). 

QLM 

(  37). 

AMASNOC  37).  CHRNO 

(  37) . REZ00420 

2 

ZP1 

(  37).  ZP2 

(  37). 

SOLID 

(  37). 

ECHCK 

(  37).  RK 

(104) .REZO0430 

3 

RL 

(  37).  RHOK 

(104). 

RDK 

(104). 

THETAK(104) .  TEMP 

(  16) .REZ00440 

4 

HEAD 

(  12).  MAXL  . 

MAXLM 

REZ00450 

c 

REZ00460 

c* 

•♦REZ00470 

C****.i******************************************************************re:Z004S0 

COMMON  /NAMEC/ 

RLF  . 

RINC 

.  RPUSH 

*  RR 

*  DR 

REZ00490 

COMMON  /NAMEC/ 

KCNT  . 

LCNT 

.  MCNT 

.  N 

.  INO  . 

IL.  IR  REZ00500 

COMMON  /NAMEC/  VAL  (  10)*  RTABRO(150>*  RTABTHI 150) *  RTABRDI 150) *REZO0510 


1  RHOTABt 150) *  THETABI150) *  RDTAB  ( 150) 
COMMON  /PROSWT/  IPROSW 
COMMON  /LMSB/  U(l)  /LMSC/  M(l) 


REZ00520 

REZ00530 


C 

DIMENSION  XZ0NE(11*9) 

DIMENSION  ZN(6) 

DIMENSION  ZK (5) *  Z(l).  WT(1) 

EQUIVALENCE  (M( 1) .KMX) *  (M(2)*Z(1))*  (M(3)»WT(1)) 
EQUIVALENCE  (RK(2) *XZONE) 

REAL  IK 
LOGICAL  1BOIL 
C 
C 

C  SET  THE  NUMBER  OF  REGIONS  TO  BE  CHANGED 

C 

NCARD  =  RK ( 1 ) 


PEZ0056J 

r<EZ00610 

REZ00600 


REZ00570 

REZ00620 

REZ00630 

REZ00640 

REZ00650 

REZ00660 

REZ00670 
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DO  140  J  a  ltNCAKO  REZ00680 

REZ00690 

INITIALIZE  VARIBLES  THAT  KEEP  TAB  ON  AREAS  TO  RE  REZONFDREZ00700 


IRN  =  0 

REZ00710 

REZ00720 

1R0  a  0 

REZ00730 

160  =  16 

REZ00740 

1L  a  XZONEdt  J) 

REZ00750 

1R  =  XZ0NE(2» J) 

REZ00760 

00  WE  TURN  ON  THE  PROFILE  SWITCH 

REZOO 770 
REZ00780 

IPROSW  a  0 

REZ00790 

REZ00800 

I F ( XR  .LT.  0)  IPROSW  a  1 

REZ00810 

1R  a  lABS(lR) 

REZ00820 

SAVRV6  a  R(IR) 

;’EZ00840 

SAVE  THE  APPROPRIATE  ARRAYS 

REZ00850 

REZ00860 

N  a  1 

REZ00870 

REZ00890 

DO  99  I  a  1L» IR 

REZ00900 

REZ00910 

REZ00920 

SAVE  TABLES 

REZ00930 

ROTAB(N)  a  RD(I) 

RTABRD(N)  a  R(I) 

REZ00940 

IF  (1.E6.1R)  60  TO  99 

RHOTAB(N)  a  l./SVd) 

REZ00970 

RE700980 

RTABRO(N)  a  R(I) 

REZ10990 

THETAB(N)  a  E(I) 

REZ01000 

REZ01010 

RTABTH(N)  a  R( I > 

REZ01020 

CAR(N)aFEWd) 

REZ01030 

CAPAC(N)aX6d) 

REZ01040 

CHIC(N)aBCd) 

REZ01050 

CHlR(N)SBRd) 

REZ01060 

X4(N)aCRTC(I) 

REZ01070 

X5(N)=RHOd) 

REZ01080 

EC(N)aGd) 

REZ01090 

IF  (EC(N).LE.O.)  60  TO  110 

PB(N)aROOd) 

REZ01100 

SMLQ ( N ) aSMLR ( I ) **2 

REZ01UO 

N  a  N  ♦  1 

REZ01120 

REZ01140 

IF  (N  ♦  1.6T.150)  60  TO  110 

.'IEZ01150 

99 

CONTINUE 

REZO 1160 

KCNT  a  IR  -  IL 

LCNT  a  KCNT 

REZ01170 

REZ01190 

MCNT  a  KCNT  ♦  1 

IF(PUSHA)  104*106*105 

REZ01210 

REZ01220 

104 

RLF  a  R0U6-1) 

REZ01230 

60  TO  111 

REZ01240 

105 

RLF  a  RD(l) 

REZO 1250 

106 

60  TO  111 

REZ01260 

REZ01270 

110 

SI  a  55.0110 

REZ01280 

CALL  UNCLE 

REZ01290 

REZ01300 

111 

1B01L  a  .FALSE. 

REZ01310 

IF(PUSHA)  112* 115*113 

REZ01320 

112 

RPUSH  a  RUB) 

REZ01330 
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IFIIL.GT.IB  .OH.  IR.LT.1B  .OH.  BOILA.EG.O.)  GO  *0  115 

REZ01340 

IBB  =  IB 

HEZ01350 

JB1  =  1BP1 

HEZ01360 

X  =  -1. 

HEZ01370 

00  TO  114 

HEZO13B0 

113 

KPUSH  -  H ( I A ) 

HEZ01390 

IF ( IL.GT. I AMI .0R.IR.LT. I AMI  .OR.  30ILA.EQ.O.)  GO  TO  115 

REZ01400 

IBB  =  I AMI 

HEZO1410 

JB1  =  I AMI 

HEZ01420 

X  3  1. 

RLZ014?0 

114 

IB0IL  =  .TRUE. 

HEZ014<40 

RPUSH=RWBl)+X*(l.-i0LID(24))*GlIBB)*S0L10I17)/A(JBl) 

HEZ01450 

115 

CONTINUE 

HEZ01460 

REZ01470 

ERROR  CHECK 

REZ01480 

REZ01490 

IFIIIR-IL)  .GT.  149)  GO  TO  100 

HEZ01500 

REZ01510 

GENERATE  VAH1-5) 

REZ01520 

REZO1530 

VAL(l)  =  0. 

REZ01540 

DO  120  I  =  2il0 

REZ01550 

VAL ( I )  s  XZ0NE ( 1+1 » J) 

REZ01560 

120 

CONTINUE 

HEZ01570 

IHM1  =  IK  -  1 

HEZ01580 

IF ( VAH t>)  .EG.  2.)  GO  TO  801 

HEZ01590 

IF(VAL(6)  .  NE.  1.)  GO  TO  100 

HEZOlbOO 

DO  800  I  =  IL* IKM1 

HEZ01610 

BOO 

VAL ( 1 )  =  VAL(l)  -  DELTARII) 

HEZ01620 

GO  TO  803 

REZ01630 

BP1 

DO  802  I  =  ILilRMl 

REZ01640 

802 

VALID  =  VAL(l)  ♦  G<1> 

REZO1650 

803 

CONTINUE 

REZO 1660 

RR  =  RIIL) 

REZ01670 

CALL  GENRAT 

REZ01680 

REZ01690 

129 

CONTINUE 

REZ01700 

HNR  =  R(NR) 

REZ01710 

RICA  =  RIICA) 

HEZ01720 

HICB  3  R(ICB) 

REZO1730 

HEZ01740 

SHIFT  ARRAYS  FOR  NEW  CONDITIONS 

REZ01750 

REZ01760 

CALL  MOVE 

REZ01770 

REZ01780 

REZ01790 

ADJUST  BOUNDARIES 

REZ01800 

IRO  =  IR 

REZ01810 

IHN  =  IL  ♦  IFIX<VAL<5>) 

HEZ01820 

IR  =  IRN 

REZ01R30 

ITEM  =  IRN-IRO 

IRM1  =  IR  1 

HI  IR)  3  SAVRV6 

REZ01840 

REZ01850 

CALCULATE  ZONE  QUANTITIES 

REZ01860 

REZ01R70 

CALL  ZONE 

HEZ01RR0 
REZO 1890 

CALCULATE  AHEA 

REZ01960 

REZ01910 

DO  1294  I  3  ILtIRMI 

REZ01920 

IF(SIO)  1291*1292*1293 

REZ01930 

1291 

All)  31. 

REZ01940 

GO  TO  1294 

HEZ01950 

1292 

AID  3  2.*  RII) 

REZ01960 

o  o  o  non 
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GO  TO  1294 

1293  All)  =  3.  *  K ( I ) **2 

1294  CONTINUE 
C 

IF  IJ.EQ.NCARD)  GO  TO  1300 

JPl  =  J  ♦  1 

00  130  I  =  JPl'NCARD 

IF  <XZONEU*J).GT.XZONE!2tI))  GO  TO  130 

XZONEdd)  =  XZ0NE(1»I )  ♦  FLOAT!  ITEM) 

XZ0NE(2.I)  =  XZ0NE(2.I)  ♦  FLOATCITEM)  *  SIGN! 1 . .XZONE! 2» I ) t 

130  CONTINUE 
1300  CONTINUE 

C 

IF(PUSHA)  136*  137*  131 

131  CONTINUE 

IB  =  IB+ITEM 
IG  =  IG+ITEM 
ICB  s  IB 

IF(IA.GT.IL)  IA  =  IS A VI ( IG» I Af RPUSH»R(1) ) 

IF(ICA.GT.IL)  ICA  =  ISAVl ( IG.-ICA.RICA.RU) ) 

NR  =  ISAVl ( IGt OfRNRtRIl ) ) 

GO  TO  139 
C 

136  CONTINUE 
IG=IG+ITEM 

IF(IB.GT.IL)  IB  =  ISAVl ( IG* -IB  *RPUSH»R ( 1 ) ) 

ICB=ISAV1 ( IG* ICBiRICUiR (1) ) 

NR  =  ISAVl ( IG»0*RNR*R(1) ) 

GO  TO  139 
C 

137  CONTINUE 

IG  =  IG+ITEM 
IB  S  IB  ♦  ITEM 

ICB  s  ISAVl ( IG*ICB»RICB»R(1) ) 

IF(ICA.GTtlL)  ICA  =  ISAVl! IG*-ICA*RICA*R! 1 J ) 

NR  =  ISAVl ( IG» 0* KNRfR ( 1 ) ) 

139  CONTINUE 
C 

I API  =  IA+1 
I AMI  =  IA-1 
IBP1  =  IU+1 
IBMl  =  IB-1 
ICAM1  =  ICA-1 
ICAP1  =  ICA+1 
ICBM1  =  ICB  -  1 
ICBP1  =  ICB  +  1 
IGM1  =  IG-1 

UPDATE  VAPOR  FRACTION  OF  BOILING  ZONE 

IF ( «N0T • IBOIL)  GO  TO  1398 
IF(PUSHA)  1392. 1398(1394 
1392  IBB  =  IB 
JB1  =  IBP1 
GO  TO  1396 
1394  IBB  =  IAM1 
JB1  s  I AMI 

1396  S0LID<24)=l.-(RPUSH-R(JBin«A(JBl)*X/(G(IBB)*S0LID(17)> 

1398  CONTINUE 


CALCULATE  NEW  VELOCITYS 

DO  200  I=IL*IR 
RDOTxR(I) 

LASTNS3 


REZ01970 

REZ01980 

REZ01990 

REZ02000 


REZ02030 

REZ02040 

REZ02050 

REZ02060 

REZ02070 

REZO2080 

REZ02090 

REZ02100 

REZ02110 

RE202120 

REZ02130 

REZ02150 

REZ02160 

REZ02170 

REZ02180 

REZ02190 

REZ02200 

REZ02220 

REZ02230 

REZ02240 

REZ02250 

REZ02260 

REZ02270 

REZ02280 

REZ02300 

REZ02310 

REZ02320 

REZ02330 

REZ02340 

REZO2350 

REZ02360 

REZ02370 

REZ02380 

REZ02390 

REZ02400 

REZ02410 

REZ02420 

REZ02430 

REZ02440 

REZO2450 

REZ02460 

REZ02470 

REZ02480 

REZ02490 

REZ02500 

REZ02510 

REZ02S20 

REZ02530 

REZ02540 

REZ02550 

REZ02560 

REZ02S70 

REZ02580 
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CALL  INTlKPCKUOT ,LASTN) 
HU<1)=RU0T 
200  C0NT1NUL 


UALANCL  TML  NLW  LNLKGIl.b 

REGENERATE  ZONE  LNEHGY,  MASS  OEPLET ION  AMU  MONFO  VARIABLES 

SM1=EC<1) 

1-IL 

bM2=G(  ID 
300  EN=0. 

1F(SUMRA,EQ,U, )  GO  TO  305 
RN=0. 

SN=0, 

305  IF(ZP1  (26)  .EQ.O, )  GO  TO  31b 
UO  310  JJ  =  1,6 
ZN(JJ)  =  0. 

310  CONTINUE 

31b  SM= AMI N1 ( SMI  * SM2 ) 

LNSEN+THETABCK) *bM 
IFlbUMRA.EQ.O.)  GO  TO  320 
RN=RN+PB(K)*SM/EC<K) 

SN=SN*SMLQ(K)*SM 
320  1F(ZP1<26).EQ.O,)GO  TO  325 
ZN<1I=ZN(1)+CAR(K)*SM 

IF  (PUSHA.GE.O,*AND,RTABRO<K) .LT.RNH)  GO  To  325 
IF  (PUSHA,LT ,  0, . AND.RTABROtK ) , GE,RNR>  GO  TO  325 
ZN(2)=ZN(2J+CAPAC<K)*SM 
ZN(3)=ZN(3)+CH1C(K)*SM 
ZN(4J=ZN(4)+CHIR(KJ*SM 
ZN(5)=ZN(5)+X4(K) *SM 
ZN(6J=ZN(6J+X5<K)*SM 
325  SM1=SM1-SM 
SM2=SM2-SM 

1F<SM1<GT<0<JGO  TO  330 
IF  (K.LT.KCNTJ  K  r  k  ♦  1 
SM1=EC(K) 

330  IF<SM2<GT<0< )GO  TO  315 
UIJrEN/GUJ 
EI(I)sEd) 

IF ( SUMRA .EQ. 0 • )  GO  TO  335 
RDO( I)=RN 

SMLR ( I )  =SQRT ( ABS ( SN/G (I)  ) ) 

335  IF(ZPl(2t>) .EO.O.)  GO  TO  340 

IF  (PUSHA,GE,0, .ANO. I.LT.NR)  GO  TO  340 
IF  ( PUSHA<LT<0< , AND, I ,GT,NR-1 )  GO  TO  340 
FEW(I)=ZN(1)/G(I) 

Xb(I)=ZN(2)/G(I) 

BC<I)=ZN<3)/G(I) 

BR<I)=ZN<4)/G<I) 

CRTC( I ) =ZN( b) /G( I ) 

RHO( I )=ZN(6) /G ( I ) 

SUMZSX6 ( I ) *WT ( 1 ) 

IF<KMX<EO< 1 )  GO  TO  337 
SUMZ=SUMZ  ♦  BC ( I ) *WT (11) 

IF (KMX>EU<2)  GO  TO  337 
SUMZ-SUMZ  +  BR( I)*WT(21) 

IF (KMX, EU. 3)  GO  TO  337 
SUMZ=SUMZ  +  CRTC(I)*WT(31) 

IF(KMX<E0<4)  GO  TO  337 
SUMZsSUMZ  *  RH05I)*WT<41) 

337  IF  ( SUMZ<EQ<0< )  GO  TO  340 

Xfe(I)  =  X6  ( I )  *  FEW(  I )  /  SUMZ 
IF  (KMX.EQ.l)  GC  TO  338 


Rt  Z02590 
HEZ02600 
.  RL/02610 
RLZO2620 
REZ02F.30 
Mt  ZO2640 
RLZ02650 
RLZO26fi0 
REZ02670 
RLZO2680 
REZ02690 
HEZ02700 
REZ02710 
REZ02720 
REZ02730 
REZO2740 
REZ02750 
REZ02760 
REZ02770 
REZC27B0 
REZO2790 
REZO2A00 
REZ 02810 
REZ02820 
REZ02830 
REZ02840 


REZ02850 

REZ02860 

REZ02870 

REZ02880 

REZ02890 

REZ02900 

REZO2910 

REZ02920 

REZO2940 

REZ02950 

REZ02960 

REZ02970 

REZ02980 

REZO2990 

REZ03000 

HEZ03010 


HEZ03020 

REZ03030 

REZ03040 

REZO3050 

REZ03060 

REZ03070 
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BC(I)  =  BC( I)  *  FEX(I)  /  SUMZ 
IF  (KMX.EQ.2)  GO  TO  338 
BR  ( I )  =  BR(I)  *  FEW(  I  )  /  SUMZ 
IF  (KMX.EQ.3)  GO  TO  338 
CRTCU)  =  CRTC(I)  *  FEW(I)  /  SUMZ 
IF  (KMX. EG. 4)  GO  TO  338 
RHO(I)  =  RHO(I)  *  FEW! I)  /  SUMZ 
•  338  CONTINUE 

340  1=1+1 
SM2=G(I) 

IF(I.LT.IR)  GO  TO  300 
CALL  UP37A ( IRO» IGO ) 

C 

C  CALCULATE  ZONE  TEMPERATURES 
C 

IF(lPROSW.EQ.l)  CALL  PROFL 
00  380  I=IL» IRM1 
IF  (ZP1 (26) .EQ.O. )  GO  TO  342 
IF  1 PUSHA • GE • 0 • • AND • I . G£ . NR )  GO  TO  364 
IF  (PUSHA. LT.O.. AND. X.LT. NR)  GO  TO  364 
342  IC  =  0 
TL=0. 

TR=2.*THETA(I“1) 

IF(I.EQ.IL)TR=2.*TH£TA(I) 

345  THET A ( I )  =  . 5* ( TL+TR ) 

SMLA(I)  =  THETA ( I ) 

CALL  EOS(I) 

IF( IC.GT.30IGO  TO  365 
OE=E(I)~EI(I) 

IF ( ABS (UE/EI ( I ) ) .LE.S19)  GO  TO  370 
IF(OE.GT.O.)  GO  TO  350 
IF(IC.GT.O)  GO  TO  355 
TL=THETA(I) 

TRSTL+2. 

GO  TO  345 
350  TR=THETA(I) 

GO  TO  360 
355  TL=THETA(I) 

360  IC=IC+1 
GO  TO  345 

365  WRITE(6» 1000)  I 

1000  FORMAT (29H0ENERGY  DISCONTINUOUS  IN  ZONE* 13) 
GO  TO  380 
364  ZK(1)  =  X6(I ) 

ZK ( 2 )  =  BC(I) 

ZK(3)  =  BR( I) 

ZK ( 4 )  =  CRTC(I) 

ZK( 5)  =  RHO(I) 

EION  =  0. 

DO  368  K  =  It  KMX 
KK  =  10*K  -  9 
LMAX  =  Z(KK) 

IK  =  AINT(ZK(K) ) 

L  =  INT(XK) 

NV  =  M(KK+3)  ♦  1 
IF  (L.EQ.O)  GO  TO  367 
00  366  JJ  =  1 »L 
NV  =  NV  +  1 

EION  =  EION  +  WT(KK)  *  U(NV) 

366  CONTINUE 

IF  (L.EG.LMAX)  GO  TO  368 

367  EION  =  EION  ♦  WT(KK)*U(NV+1)+(ZK(K)-IK) 

368  CONTINUE 

369  AMASNI  =  AMASN(I) 

CV(I)  =  1.5  •  AMASNI  *  (1.0  +  FEW(D) 


REZ03120 

REZ03130 

REZ03140 

REZ03460 

REZ03150 

REZ03160 

REZ03170 

REZ03180 

REZ03190 


REZ03210 

REZ03220 

REZ03230 

REZ03240 

REZ03250 

REZ03260 

REZ03270 

REZ03280 

REZ03290 

REZ03300 

REZ03310 

REZ03320 

REZ03330 

REZ03340 

REZ03350 

REZ03360 

REZ03370 

REZ03380 

REZ03390 

REZ03400 

REZ034.0 
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IHETA(I)  =  (Ed)  -  AMASNI*EI0.M  -  EZERO(I>>  /  CV(I>  02/23/66 

P1(I>  =  AMASNI  *  (1.0  *  FEW( l ) )  *  THETA ( I )  /  SV(I) 

370  EI(I)=£(I)  REZ03420 

380  CONTINUE  REZ03430 

REZ03440 
REZ03450 
REZO3470 

140  CONTINUE  REZO34B0 

RETURN  REZO3490 

C  REZ03500 

C  ERROR  REZ03510 

100  SI  S  55.0100  REZ03520 

CAUL  UNCLE  •  REZ03530 

0  REZ03540 

END  REZ03550 
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SUUROUIINE  CONVIN 
C**  t  COMPILED  JUNE  (>#  I9w> 


C* 

C* 

c* 


SPUTTER  C  0  M  M  0  N 


COMMON  LMDA(940b) 

DIMENSION  EC(lb2) 

EQUIVALENCE  (SI#  LMOA(lJb))#  (ECU),  LMDA(SM7) ) 


COPVOOIO 

♦♦♦♦♦convorvo 
♦copvno  v> 

♦♦COPV0040 

♦  *coi)vnn',o 

♦  ♦convnofio 


c* 


v.  vi  i  luunv 

ClltMtMMMMIMMttMtMMtmHttitttm*..#*.*..  .  ^  **CONV0090 

COMMON  /NAMEC/  MC#  .  ,«NC  .»#»,.»„  » . SRS!! 

KCNT  >  LCNT  •  MCNT  •  n  •  IMO  #  It#  tr  rnMunl-in 

THETABCiSt!;  StS°;KSJ#  RWoasoi.CJlVOjSo 

I/SWT  CONVOlbO 


COMMON  /NAMEC/ 
COMMON  /NAMEC/ 
i  RHOTABU50)* 


COMMON  /SREZS/ 

DOUBLE  PRECISION  R0#R1 #R2#HJ,G1 #G2.VS#/M!. 


M 

X 

X 

X 

0 

0 

0 

X 


Ml  MN  ( SPACE  COORDINATES  GIVEN) 

X  X  SPECIAL  CASE  (X  MEANS  VALUE  GIVIN) 

0  X  SPECIAL  CASL 

X  X  ERROR 

o  X  ehror 

X  0 

0  0 


c 


MASS  -  SPACE 
IV6  =  VAL (6)  ♦  0.1 


SWITCH 


C 

C 

C 

C 

C 

C 

C 

C 

C 


ITEM  =  0 

TEM  =  ABS(VAL(1) ) -Al)S(  VAL  ( 3) J 
DO  190  I  =  1,3 


IF<  (VAL(b)-l,5)*VAL(I)  .GE.O. )  GO  TO  190 
RO=DBLE(RR) 

K=MIN0(KCNT#2) 

80  IF(RTAURO(K) .GT.RR.OK.K.EG.KCNT)  GO  TO  90 
k=k+i 

GO  TO  80 
90  CONTINUE 

GO  TO  (110,150) • IVb 


SPACE  COORDINATES 


110  CONTINUE 


CHECK  FOR  MIS-SPECIFIC AT ION 

ERROR 

IF(VAL(3) (GT.0.0  .AND.  VAL(l) .EQ.O.O)  GO  TO  180 
CHECK  FOR  SPECIAL  CASE 


120  CONTINUE^  ’AND*  (VALa,‘GT*°-°  -AND.  VAL (3) .GT.0.0) )  GO  TO  130 
RlsRO 

IF(KCNT.GT.2)  GO  TO  122 
G1  =  1.D20 
GO  TO  124 

122  !ZS"T-E0-  '-*  * »  «  « *« 

G1=ZMS(R1,R2,k) 

GO  TO  124 

123  GlsOOLE(EC(K-l) ) 

124  G2SDULE ( VAL ( I ) ) 

IFIG2.GE.G1 )  GO  TO  126 


C0NV0180 

CONV0190 

COMV0200 

CONV0210 

CONV0220 

CONV0230 

CONV0240 

CONV0250 

CONV0260 

CONV0270 

C0NV0280 

C0NV0290 

CONV0300 

CCNV0310 

C0NVC320 

CONV0330 


CONV0350 

CONV0360 

CONV0370 

CONV0410 

CONV0420 

CONV0430 

CONV0440 

COMV0450 

CONV0460 

CONV0470 

CONV0480 

CONV0490 

CONV0500 
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c 


c 


c 


c 


c 


12b  CAL-  ZKALHRl#  K3»  02 »  V j. K ) 

00  TO  128 

12b  IF(G2.GT.G1>  GO  TO  127 
R3=R2 
00  TO  128 

127  R1=R2 
02=G2-G1 

IFIK.EQ.KCNT)  GO  TO  125 

K=K+1 

GO  TO  122 

128  VAL( I)=SNGL(R3-R0) 

GO  TO  190 

130  CONTINUt 

VAL(3)  =  TEM 
ITEM  =  1 
GO  TO  120 

MASS  COORDINATES 

150  CONTINUt 

CHECK  FOR  MIS-SPECIFIC ATI ON 

ERROR 

IF(VAL(3) .LT.0.0  .AND.  VAL( 1) .EG. 0 . 0 )  GO  TO  lftO 
CHECK  FOR  SPECIAL  CASE 

IF ( I .EQ. 3  .AND.  <VAL<1>. LT.0.0  .AND.  VAL13) .LT.0.0) )  GO  TO  17r 

160  G1  =  0.00 
Ro=R0+UbLE  ( ABS  ( V AL  ( 1 1)  1 
R1=R0 

IF1KCNT .GE.2'  GO  TO  162 

161  G1=ZMS(R1»R3.K)+Gl 
GO  TO  165 

162  R2=0bLE(RTAbR0(K>) 

IF(R2.GT.K3)  GO  TO  161 

IF  (SNGLtRl )  .  EQ  .RTAURO  (K-l  1 .  AND.  IZSWT  .EO.(-D)  GO  TO  163 

G1=ZMS(R1.R2.K)+G1 

GO  TO  164 

163  Gl=DbLE(EC<K-l))+Gl 

164  IF(R2.LQ.R3>  GO  TO  165 
R1=R2 

IF (K .EU.KCNT)  GO  TO  161 

K=KU 

GO  TO  162 

165  VAL(I)  =  -SNGL(Gl) 

GO  TO  190 

170  continue 

V’L(I)  =  TEM 
item  =  1 

GO  TO  160 

180  continue 

181  SI  =  43.0181 
WRITE  <6»1P) 

RETURN 

190  CONTINUE 

IF ( ITEM  .EQ.  0)  RETURN 

VAL(3)  =  AbS(VAL(l) )  -  ABS(VALC3>) 

RETURN 


CONV0590 

CONV0600 

CONV0610 

CONV0620 

CONV0630 

CONV0640 

CONV0650 

CONV0660 

CONV0710 

CONV0720 

CONV0730 

COMV0740 

COMV0750 

CONV0760 

CONV0770 

CONV0780 

CONV0790 

CONVOBOO 


CONV0870 

CONV0880 

CONV0890 

CONV0900 

CONV0910 

CONV0920 

CONV0930 

CONV0940 

COHV0950 

CONV0960 

CONV0970 

_ CONV09B0 

CONV0990 

CONV1000 

CONV1010 

CONV1020 


10  FORMAT  (  52H1A  MIS-SPECIFICATION  HAS  HEEN  FOUND  IN  TIIF.  ZONE  DATA)  CONV1030 
END  CONV1040 
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C  COMPlLtD  JULY  15.  t'i(ni 

C4****i*  ♦♦♦***♦$1**44*444  4444  *4*4*44  44  14  44  4444  44*44*44  444  |  4444  44444444444 

i>  H  U  T  1  |  l<  (.  O  M  '1  <1  ||  ,! 


COMMON 

1  KMAX 

2  II 

3  IIP1 

4  FREU 

5  CVA 

b  EMINA 
7  RHOL 
B  RPI A 
COMMON 

1  DTMAX1 

2  GAMA 

3  ROIA 
F  S4 

!>  S12 

i  S20 
7  ASO  ( 
)  DELR  ( 


LMDA ( 37 )  #  Nl< 


BLANK  1 , 
IG  , 
IGM1  , 


CNTMAX#  AH 


BLANKS* 
NKAD  * 
I ALPHA* 


CVU 

EMINU 

RHOR 

RPIB 

TD 


SLUG 

CA 

LPIO 

RPDIA 

TE 


NSMI  It  .  I A  .  I|t 

III. AUK  3.  I  API  .  IHPI 

ULAIIK4  *  |AM!  *  Iltf'I 

uLAtlK!)*  Tli  *  Imax 

ASMI.R  *  PUSIIA  *  PI  if .hu 

ALPHA  .  HVA  ,  MVP 

CH  *  GA  ,  (,M 

LPSI  *  RIA  ,  KIH 

PPDIII  * 

DTII?  * 


UTMAX2*  DTMAX3.  UTK 


WCRIT  *  SIGMAGi 
ROIAM1*  Roll)  > 


COMMON 


*  S5 

*  S13 

*  EO 
(152)* 
(  37)* 


(152)*  RU 
(  37)*  P 
P2 

W  (152)*  E 
V  (152)*  G 
X3  (152)*  X4 
SMLA  (152).  SMLH 
EC  (152)*  EH 
BIGB  (152)*  CV 
CHIR  (152).  CAPAC 
CRTPC  (152)*  GOFR 


AC 

K01HP1 

57 

Sib  1 
TAU  , 


Sb  ,  S7 

S14  ,  Sib 

FO  .  TAU 

(lb2).  vn 

(152).  PI 
( lb2 ) *  SV 
(152)*  El 
(152)*  0 
(152)*  X5 
(152).  SMLC 
(152)*  SMLO 
(152).  HC 


TPHIIJT*  TA 
OTH2P  *  l)T|l] 
SWITCH*  CO 
AC03T4 *  CNVRT 
GMS  *  SI 


*  SO  *  S9 

*  S16  »  SI" 

*  ZERO  *  H 
(152)  *  ROD 
(152)*  PH 
(152)*  RHO 
(152)*  1.K 
(152).  C 
(152)*  Xft 
(152)*  SMLt) 
(152).  SMLH 
(152).  MR 


1(  A 
I(  API 
H  AMI 
I<LAIII<<1 
ItOIl  A 
MCA 
GL 

RDIA 

Til 

OTRMIN 

CHIN 

SUMRA 

S2 

sin 

sin 


*  * 
* « 
♦  ♦ 

I  CM  * 
I CtIP l  * 
1CIIW1  , 
lil.LPHT  * 
non  m  , 

HCIj  * 

gu  , 
koim  * 

TC 

OTMAX  * 
DELTA  * 
SUMP II  * 
S3  * 
Sll  * 
519  , 


(152),  DELTAR  ( lb2)  • 
(152),  SKLH  (152), 
(152),  rill  (152) 
(152),  THETA  (152), 


(152).  CAI'AR  (152),  CRTC 


COMMON 

1  FRLM 

2  ZP1 

3  RL 

4  HEAD 


(  37),  ZP2 
(  37),  RHOK 
(  12),  MAXI 


(152).  FEW 
(37).  EKLM 
(  37).  QLM 


(152),  CAR 
(  37),  ELM’ 


f  1  52 ) ,  A 
(152),  X2 
(152),  X7 
(152),  SMLE 
0  52)5  MIGA 
(152),  C)I1C 
(152),  CRTR 
(152),  OKLM 
(  37),  FCLM 


(  37),  SOLID  (  37), 
(104),  ROK  (104), 
*  MAXLM 


A  (152), 
X2  (152), 
X7  (152), 
SMLE  (152). 
RIGA  052), 
C)IIC  (152), 
CRTR  (152), 
OKLM  (  37) 
FCLM  (  37), 
CHHNO  (  37), 


(  37),  AMASNOI  37),  CHHNO  (  37), 
(  37),  ECMCK  (  37),  RK  OOP), 
(104),  THLTAK (104),  TEMP  {  16), 


. . . . * . 


COMMON  /NAMEC/ 
COMMON  /NAMEC/ 
COMMON  /NAMEC/  VAL 


RLF  ,  RINC  ,  RPUSH  ,  RR  ,  DR 

KCNT  ,  LCNT  ,  MCNT  ,  N  ,  INO  *  IL,  IR 

(  10),  HTAURO( 150) ,  RTAHTIK  150) ,  RTAUROI 150) , 


1  RHOTAB(150),  THETAB(lbO),  RDTAU  (150) 

Ct  JMON  /SREZS/  IZSWT 

i.  <vwLE  PRECISION  R0*H1*H2,H3,R4,G1,G2 ,GS, VS, VOL , 7MS 
CALCULATE  G( I / ,R( I ) , DELTAR ( I ) , SV ( I ) 

IRM1SIR-1 
IRM2.=  IR  -  2 
ILP1=IL+1 
K=MIN0(KCNT,2) 

80  IF(RTABRO(K) ,GT,R( IL) .OR.K.EW.KCNT)  GO  TC  90 
K  — K  ■*  1 
GO  TO  80 
90  CONTINUE 

IF(VAL(6) .EQ.2.)  GO  TO  120 
SPACE  ZONING 
CALCULATE  ZONE  WIDTHS 
DELTAR (IL)SVAL( 2) 

IF  ( ILP1.GE, IRM2)  GO  TO  930 
DO  92  I  =  ILP1.IRM2 
t)ELTAR(I)=DELTAR(I-l)*VAL(4) 

92  CONTINUE 

DELTAR ( IRM1)  s  VAL(3) 

DO  93  I  =  ILP1.IRM2 
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J  =  ILP1  ♦  IRM2  -  I 

WTR  =  FLOAT ( IHM1  -  J)/H  OAT< IRM1  -  »L) 

UELTAR(J)  =  WTR+ULLTAIM  J)  ♦  ( 1  .-WTR)  *UI.LI  AIM  J+ 1  )/VAl  (M) 

93  CONTINUE 

C  CALCULATE  ZONt  BOUNDARIES 

SOU  Ml  =  0BLE(R( IL) ) 

DO  94  1  :  ILPl*Ii<Ml 

Ml  =  HI  +  DBLE (ULLTAK ( 1-1 ) ) 

1M I )  =  SNGL(Rl) 

94  CONTINUE 

C  INITIALIZE  LEFT  HOUNUAKY  OF  ZONING  REGION 

I  =  IL 

KQ=DULE(R( IL) ) 

104  K4rDBLE(HTABR0(K) ) 

C  INITIALIZE  ZONE  I 

105  GS  =  0.D0 
RlsRO 

H3=R0+DbLE (DELTAR (  I) ) 

C  SUBZONE  CALCULATION 

106  R2=R3 

IF(K.NE.KCNT)  H2=0MIN1 (H3*R4) 

GS=ZMS(R1*R2*K)  +  GS 
K1=R2 

IFIR2.LT.R4. OR. K.LQ.KCNT)  GO  TO  100 
C  INITIALIZE  NEW  SECTION  OF  DENSITY  PROFILE 
K=K+1 

R4=DULE (RTABROIK ) ) 

C  IS  ZONE  I  COMPLETE 
108  IFIR2.LT. R3)  GO  TO  106 

C  ESTABLISH  MASS*  SPECIFIC  VOLUME  AND  RIGHT  BOUNDARY  OF  ZONE 
G'I)sSNGL(GS) 

VOL=R3-R0 

GO  TO  (114*111*112)* 1 ALPHA 

111  VOL=VOL*(R3+RO) 

GO  TO  114 

112  VOL=VOL* (R3+*2+R3*R0+R0**2) 

114  SV ( I ) =SNGL ( VOL/GS ) 

1  =  1  +  1 

IF  (I.EQ.IR)  GO  TO  131 
C  GENERATE  NEW  ZONE 

R0=R3 
GO  TO  105 
C  MASS  ZONING 

120  CONTINUE 

C  CALCULATE  ZONE  MASSES 

G( IL)=VAL(2) 

DO  121  I  =  ILP1 * IRM2 
G ( I ) =G ( 1-1 ) +VAL (4  > 

121  CONTINUE 
G(IRMl)  =  VALI3) 

DO  122  I  =  ILP1 »  IRM2 
J  =  ILP1  +  IRM2  -  I 

WTR  =  FLOAT ( IRM1  -  J)/FL0AT(IRM1  -  IL) 

G(J)  ?  WTR*G( J)  +  ( 1 • -WTR ) +G ( J+l ) /VAL ( 4 ) 

122  CONTINUE 

C  INITIALIZE  LEFT  BOUNDARY  OF  ZONING  REGION 

K3=DBLE(H( IL) ) 
plL 

IF  (K.LT.KCNT)  GO  TO  123 
G1  =  1.D20 
GO  TO  125 

123  IF(R(IL) .EQ.RTABRO(K-l) .AND.IZSWT.EQ.  (-1) )  GO  TO  124 
Gl=ZMS(R3*DBLE(RTABRO(K) ) *K) 

GO  TO  125 

124  G1  =  DBLE  (EC  (K"*!) ) 

C  INITIALIZE  ZONE  I 

125  G2=DBLE(G( I) ) 

R2=R3 
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vol  =  0.00 

C  SUBZONE  CALCULATION 
12b  GS=DMIN1 (G1 »G2) 

K1=R2 

CALL  ZRJ«D(Rl»R2i6St  VStK) 

GlsGl-GS 

G2=82-GS 

VOL=VOL+VS 

IF  (G1.6T.0.D0)  GO  TO  127 
IF  (K.LT.KCNT)  GO  TO  1260 
G1  =  G2 
GO  TO  127 

C  INITIALIZE  NEK  SECTION  OF  DENSITY  P..OFILL 
1260  K  =  K  +  1 

G1=DBLE(EC(K-1)> 

IF ( IZSKT.EQ.O )  Gl=ZMS(DBLE(RTAHKO(K-l) S  «DUl  L(RTAnRO(K) )  >K) 
C  IS  ZONE  I  COMPLETE 

127  IF  (G2.GT.0.D0)  GO  TO  126 

C  ESTABLISH  WIDTH*  SPECIFIC  VOLUME  AND  RIGHT  BOUNOARY  OF  ZON) 
DEL TAR ( I ) -SNGL ( R2-R3 ) 

R3*R2 

SV ( I ) =ONGL ( VOL ) /G ( 1 ) 

1=1+1 

IF  (I.EQ.IR)  GO  TO  130 
R(  I  )="*NGL<R3) 

GO  TO  125 

C  IS  RIGHT  REGION  BOUNDARY  PRESCRIBED 

130  IF  (R(IR).£Q.1.E20)  GO  TO  ISO 

131  R2  s  DBLEiR(IR) ) 

DO  135  I  s  ILP1 » IRM1 

J  =  ILP1  ♦  IRM1  -  l 

KTR  r  FLOAT ( IR-J ) /FLOAT ( IR“IL) 

R2  =  R2  -  DBLECDELTAR(J)) 

R(d)  =  MTR+RIJ)  +  ( 1 • "KTR ) +SNGL ( R2 ) 

135  CONTINUE 
GO  TO  160 

150  R(IR)  =  SNGLIR3) 

160  RETURN 
END 
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SUBROUTINE  ZRA0(X1*X2*GG* VV*KK> 

COMMON  LMOA(9405) 

EQUIVALENCE  « SI  *  LMO  A  (1 36) )  .  it  ALPHA * LMDA  ( 82 ) ) 
C************* **************************************** 
COMMON  /NAMEC/  RLE  i  RINC  »  RPUSH  ♦  RH 

COMMON  /NAMEC/  KCNT  t  LCNT  #  MCNT  *  N 

COMMON  /NAMEC/  VAL  (  10) «  RTABRO(150>*  RTABTHC 
1  RHOTAB(150)  »  THETAB(  lF'.'i  #  RDTAB  ( 150) 

C 

C  COMPUTE  RIGHT  BOUNDARY  R2  ANO -VOLUME  VS  GIVEN  LE 
C  ANO  MASS  GS  OF  A  REGION  IN  WHICH  DENSITY  IS  LINE 
COMMON  /SREZS/IZSWT 

OOUBLE  PRECISION  R1»R2»GS#VS#C1#C2»X1»X2»GG#VV 

DOUBLE  PRECISION  C3»  C4*  Cb#  Co 

RlaXl 

GS-GG 

KsKK 

IFtK.GE.2)  G  J  TO  5 
K=2 

GO  TO  12 

5  IF  (RTABRO(K).NE.RTABRO(K-l))  GO  TO  10 
R2-R1 
VS  =  O.DU 
GO  TO  60 

10  IF  ( IZSWT .EG. (-1 ) )  GO  TO  12 
S  s  RHOTAB(K)  -  RHOTAB(K-l) 

IF(S.NE.O.)GO  TO  3C 
12  VS  S  GS/DBLE(RHOTAB(K-D) 

GO  TO  (15* 20 i 25 )» I ALPHA 
15  R2  s  VS  ♦  Rl 
GO  TO  60 

20  R2  S  OSQRT(VS  *  Rl**2 ) 

GO  TO  60 

2b  R2  a  (VS  ♦  R1**3)**(1.D0/3.D0) 

GO  TO  60 

30  A1  a  RTABRO(K-l) 

A2  a  RTABRO(K) 

Cl  =  Rl 

GO  TO  (34*32*33)*  I ALPHA 

32  A1  a  Al**2 
A2  a  A2**2 
Cl  a  Rl**2 
GO  TO 

33  A1  a  Al**3 
A2  a  A2**3 
Cl  a  Rl**3 

34  S  a  S/(A2  -  Al) 

C3  a  DBLE(A1  -  RHOTAB (K-l ) /S) 

C4  a  GS*DBLE (2./S) 

Cb  a  Cl  -  C3 
C2  a  C3 

C6  a  C5**2  +  C4 

IF  (C6.GTt0.00)  C2  a  C3  ♦  DSIGN(1.DO*C5)*DSORT(C 
GO  TO  (35*36*37)*  I ALPHA 

35  R2  a  C2 
GO  TO  38 

36  R2  a  DSQRT ( C 2 ) 

GO  TO  38 

37  R2  a  C2**(l •00/3*00) 

38  VS  a  C2  -  Cl 
60  X2SR2 

WaVS 

RETURN 

END 


DR 

INO  *  IL*  IR 
0)*  RTARRD(ISO)* 


BOUNDARY  Rl 
IN  R**ALPHA 
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00UF»LE  PRECISION  FUNCTION  2MS(R1  *R2»KK) 

C  COMPILED  APRIL  lb'  1966 

C*** ****************** *********** ************* ************************** 

COMMON  /NAMEC/  RLF  »  HINC  *  KPUSH  *  HR  *  UR 

COMMON  /NAMEC/  KCNT  »  LCNT  *  MCNT  *  N  *  INO  ,  IL'  IP 

COMMON  /NAMEC/  VAL  <  10) »  RTAHRO(lbO) »  HTAI1TH( ISO) t  RTAHRIX  ISO) » 
1  RHOTAU< ISO) *  THETAH(lbO)*  RUT Al)  (ISO) 

FORMERLY  NAMEU  SAV 

COMPUTES  MASS  FOR  A  REGION  WITH  DENSITY  LINLAR  <IZSWT=0>  OR 
CONSTANT  <IZSWT=-1)  WITH  RESPECT  TO  R** I  ALPHA 
COMMON  LMOA(9405)  N 

EOUl VALENCE (SI 'LMUA(lJb) ) » ( I  ALPHA »LMDA( 62) )  * ( ALPHA #LMDA( 79) ) 

DOUBLE  PRECISION  Rl»R2»ZMS*Ul *H2*R1A»H2A 
COMMON/SREZS/IZSWT 
K=KK 

IF(K.GE.2)  GO  TO  b 
K-2 
2  S=0. 

GO  TO  12 
5  ZMS  =  0.U0 

IF  (RTAURO(K) .EQ.HTAURO(K-l) )  RETURN 
IF  (IZSWT.EQ.(-li)  GO  TO  2 
S  =  RHOTAB(K)  -  HHOTAlXK-1) 

12  R1A  S  R1 
R2A  =  R2 
A1  =  RTABRO (K-l ) 

A2  =  RTABRO(K) 

GO  TO  (25' 15*20) ' 1  ALPHA 
lb  H1A=R1A**2 
R2A=R2A**2 
A1  =  Al**2 
A2  =  A2**2 
GO  TO  2b 
20  R1A=R1A**3 
R2A=R2A**3 
A1  =  Al**3 
A2  s  A2**3 
2b  S  s  S/(A2-A1) 

B1  i  UBLE <  SHOT AB ( K-l )  -  S*A1) 

B2  s  DBLE(.5  *  S) 

ZMS  s  Ul*(r<2A-RlA)  ♦  B2*(R2A**2  -  R1A**2> 

IF(ZMS.GE.O.OO)  RETURN 

SU48.0025 

CALL  UNCLE 

END 
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SUBROUTINE  PROI-L 

. * . . 

UTTER  COMMON  ** 


C 

C* 

C* 

C* 


S  P 


COMMON 
KMAX 


LMDA(37). 
ULANK1. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


2 

11 

9 

IG  * 

3 

IIP1 

9 

IGM1  * 

4 

FREQ 

9 

CNTMAX. 

5 

CVA 

9 

CVB  * 

6 

LMINA 

9 

EMINB  . 

7 

RHOL 

9 

RHOR  » 

8 

RPIA 

9 

RPIB  * 

COMMON 

TO  • 

1 

DTMAX 1. 

0TMAX2* 

2 

GAMA 

WCRIT  . 

3 

ROI A 

ROI AMI . 

4 

54 

S5  » 

5 

S12 

S13  * 

6 

S20 

EO  * 

7 

ASQ 

(152).  RD 

8 

DELR 

( 

37).  P 

COMMON 
1  W 
V 

X3 

SMLA 
EC 

BIGB 
CHIR 
CRTPC 
COMMON 

1  FRLM 

2  ZP1 

3  RL 

4  HEAD 


NR  .  NSMLR 
BLANK2.  BLANKS 
NR AD  .  BLANK4 
I ALPHA i  BLANKS 
ASMLR 
ALPHA 
CB 

EPSI 
RPDIB 
DTH2 


AR 

SLUG 
CA 

EPIO 
RPDIA 

T 

D.MAX3.  DTR 
SIGMAQ.  AC 


ROIB  » 
S6  * 

S14  * 

FO  » 

(152)* 


KOIBP1 

S7 

SIS 

TAU 

VD 


1A 

I  API 
IAM1 
TH 

PUSHA 
HVA 
GA  ’ 
RIA  * 
TPRINT * 
DTH2P  . 
SWITCH* 
AC03T4. 


IB 

IBP1 

IBM1 

TMAX 

PUSHB 

HVB 

GB 

RIB 

TA 

DTH1 

CO 

CNVRT 
SI 
S9 
S17 
R 


152)i 
152)* 
152)  • 
152)* 
152)* 
152)* 
152). 
152). 

37). 

37). 

37). 

12). 


E 

G 

X4 

SMLB 
ER 
C  V 

CAPAC 

GOFR 

TELM 

WLM 

ZP2 

RHOK 

MAXL 


(152).  PI 
(152).  SV 
(152).  El 
(152).  D 
(152).  X5 
(152).  SMLC 
(152).  SMLQ 
(152).  BC 


GMS 

se 
S16 
ZERO 
(152). 

(152). 

(152). 

(152). 

(152). 

(152). 

(152).  SMLD 
(152).  SMLH 
(152). 


ICA  . 
ICAP1  . 
ICAMl  . 
BLANKb. 
BOILA  . 
HCA  . 
GL  * 
RDIA  . 
TB  * 
DTRMIN. 
CMIN  . 
SUMRA  . 
S2  . 
S10  * 
SIB 


* 

** 
** 
** 

ICB  . 
ICBP1  . 
ICRM1  . 

delprt. 

B01LB 

HCB 

GR 

RDIB 

TC 

DTMAX 
DELTA 
SUMRB 
S3 
Sll 
S19 


RDD 

PB 

RHO 

EK 

C 

X6 


BR 


l  □  a*  — 

(152).  CAPAR  (152).  CRTC 

(152).  FEW  (152).  CAR 

(  37).  ELM 


(152).  DELTAR( 152) * 
(152).  SMLR  (152). 
(152).  PB1  (152) 
(152).  THETA  (152). 
(152).  A  (152). 
(152).  X2  (152). 
(152).  X7  (152). 
(152).  SMLE  (152). 
(152).  RIGA  (152). 
(152).  CHIC  (152). 
(152).  CRTR  (152). 
(152).  OKLM  (  37) 

(  37).  FCLM  (  37). 


(  37 1 9  fcL.ro  »  . 

<  3T> *  *MASNO(  37).  CHRNO  37  . 


(37).  EKLM 

1  I]]]  5SlO  »  3T»;  ECHCK* (  37).  RK 
(104).  ROK  (104).  THETAK( 104) .  TEMP 
.  MAXLM 


(104). 
(  1G)» 


•* 


COMMON  /NAMEC/ 
COMMON  /NAMEC/ 
COMMON  /NAMEC/  VAL 

_ _  .  .  T  l  lC 


KLF  *  RINC  .  RPUSH*  “ 

KCNT  •  LCNT  .  MCNT 
(  JO).  RTABROI 150) . 


.  RR  *  DR  ..  „ 

,  N  .  INO  .  IL.  IP 

RTABTH(150).  RTABRD( 150) . 


i  ;CeT*eiijo>.  aoT*B  uso. 

SOLVE  QUADRATIC  FOR  COEFFICIENTS 

El  s  AO  ♦  A1*R1  ♦  A2»R1**2 
E2  s  AO  ♦  A1*R2  ♦  A2*R2**2 


ZSKR2.R1) 

ZS2(R2.R1) 

ZS3(R2*R1> 

ZS4(R2.R1) 

ZS5(R2.R1» 


=  R2-R1 

s  R2  ♦  R1  .  „ 

-  R2**2  ♦  R1*R2  ♦  Rl**2 

l  *  **- 


IRMl  =  IR-1 
EG  =  0.0 
K  =  IL 
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00  100  I  =  IL» IHMl 
IF  ( £ ( I ) . NE • 0 1 )  00  TO  20 
K  =  I  ♦  1 
GO  TO  100 

20  EG  =  EG  ♦  E(I)*G(I) 

100  CONTINUE 

IF  IK.GE.IR-2)  RETURN 
El  =  E(K-l) 

E2  ~  EUR) 

HI  =  R (K) 

R2  =  R(IR) 

10'  CONTINUE 

EGS  =  SV (K)  *  EG 

CALCULATE  COEFFICIENTS 

U1  =  ZS1 (H2»R1) 

GO  TO  «110»120#130)#IALPHA 
110  CONTINUE 

B2  =  (0.5*ZS2(R2'R1)-R2) 

B3  =  (ZS3(R2»R1 )/3.0-R2**2) 

GO  J0(}5SS/B1“E2,*<E2“E1,*B2/B1,/I_ZS2,R2'RH*B2+B3) 

120  CONTINUE 

B2  S  (2,0*ZS3(R2»R1 )/ (3. 0*ZS2(R2»R1 ) )-R2) 

B3  s  0.5»ZS4(R2.R1)/ZS2IR2*R1)-R2**2 

GO  f0(J5SS/<B1*ZS2(R2’Rln‘E2,"<E2“E1,*B2/B1,-'«-2S2(R2»Rl)*B2*R3) 
130  CONTINUE 

B2  =  <0.75*ZS4<R2,R1)/ZS3(R2,R1)-R2) 

B3  =  <0.60*ZS5(R2.R1)/ZS4«R2,R1)-R2**2) 

A2  =  I < EGS/ ( B1 *ZS3 ( R2 . R  J ) ) -E2 ) - 1 E2-E1 ) *B2/B1 ) / ( -ZS2 (R2»R1 ) *B2*83) 
140  CONTINUE 

A1  =  ( E2“E1 ) /B1 -A2*ZS2 ( R2 » R 1 ) 

AO  s  E1-A1*R1-A2*R1**2 


FIND  THE  ENERGY  FOR  EACH  ZONE 


DO  190  I  =  K#  IRM1 
R2  =  RU+l) 

HI  =  R ( I ) 

C 

GO  TO  ( 150* 160*170) * IALPHA 
ISO  CONTINUE 
C 


C 


E«D  =  A0+0.5*A1*ZS2IR2.RU  ♦  A2*ZS3(R2»Rl)/3.0 
GO  TO  180 


c 


160  CONTINUE 

*  2»0*Al*ZS3(R2*Rl)/3,0  ♦  A2*ZS4(R2,Rl)/2.0 

GO  TO  ISO 


170  CONTINUE 

Ed)  =  AO  ♦  0  •  75  *  A  1*ZS4  ( R2 » R 1 ) 
180  CONTINUE 
EKI)  s  Ell) 

ICU  =  I 

IF  (EU).GE.El)  SO  TO  190 
Gp  TO  199 
190  CONTINUE 
RETURN 


♦  0 .60*A2*ZS5(R2»R1 ) 
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199  CONTINUE 

DO  200  J=K»ICU 

UJlsEl 

El(J)s£(J) 

LG=EG-E( J) *G(J) 

200  CONTINUE 
202  K  s  ICU  ♦  1 

IFCCU.EQ.IRH1)  RETURN 
HI  =  R<K) 

K2  =  R(IR> 

GO  TO  101 

ENO 
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SECTION  II 

FORTRAN  ANALYSIS  ROUTINE 
2.1.  INTRODUCTION 

The  FORTRAN  Analysis  Routine  i.  written  entirely  in  FORTRAN  IV 
for  a  65K-UNIVAC  1108  computer.  Conversion  to  other  dialects  of 
FORTRAN  IV  is  not  expected  to  present  major  obstacles.  This  routine 
constats  of  two  independent  programs.  Program  One  is  a  variable 
renaming  program;  i.e.,  it  provides  for  replacement  of  specified  variable 
name,  by  other  names  in  a  FORTRAN  source  deck,  edit,  ,h.  deck  by  pro. 

g  ppropnate  spacing,  punches  out  a  new  sequenced  deck,  and  provides 
a  t sting  of  this  deck.  In  addition,  it  generate,  a  tap.  containing  ali  the 
required  information  for  the  other  program.  Program  Two  provides  a 
printout  of  an  extensive  cross  reference  dictionary  for  the  variable,  used 
.0  th.  program  and  the  source  line  in  which  each  variable  appeared,  in- 
eluding  a  description  of  the  variable,  it.  use,  in  each  subprogram,  and  the 

subprogram  in  which  it  appears.  It  is  intended  that  these  programs  be 
widely  used  for  program  maintenance. 

2.2.  DESCRIPTION 

M  general,  the  input  to  Program  On.  consists  of  an  unlimited  number 
of  compilable  FORTRAN  IV  subprograms,  each  of  which  is  preceded  by  up 
to  15  "valid"  string,  which  may  appear  in  the  subprogram,  that  follow. 

For  each  "valid"  string  to  be  substituted  for,  a  substitution  string  must  be 
provided  which  maintains  the  compilability  of  th,  newly  created  statement. 

A  String  is  a  sequence  of  characters.  Only  a  valid  string  may  be 
modified.  To  be  valid,  a  string  must  be  one  of  the  following. 
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1.  A  variable  name  or  a  common  block  name  or  a  subprogram 
name  (hereafter  the  word  "identifier"  will  be  used  to  refer 
to  one  of  these  three), 

2.  A  variable  or  subprogram  name  followed  by  a  balanced  set  of 
parentheses, 

3.  The  contents  of  a  balanced  set  of  parentheses, 

4.  Any  string  contained  between  consecutive  commas, 

5.  Any  string  contained  between  a  left  parenthesis  and  the  first 
comma  that  follows  it,  if  any,  prior  to  the  corresponding  right 
parenthesis, 

S' 

6.  Any  string  contained  between  che  right-most  comma  following 
a  left  parenthesis,  if  any,  and  the  corresponding  right  paren- 
the  sis; 

and  it  must  not  be  any  of  the  following: 

1.  An  arithmetic  or  logical  operator, 

2.  A  single  non-alphabetic  character, 

3.  .TRUE.,  .FALSE.,  .T.,  .F.  , 

4.  A  real  or  double  precision  or  integer  constant, 

5.  Any  hollerith  string  (as  appearing  in  a  format  statement  or 

data  statement  or  as  an  argument  of  a  call).  ^ 

In  an  arithmetic  statement  or  arithmetic  function  statement,  the  left 
side  of  an  equal  sign  is  always  a  valid  string;  however,  the  right  side  may 
or  may  not  be  a  valid  string.  The  entire  contents  of  the  balanced  set  of 
parentheses  following  an  IF  statement  is  not  necessarily  a  valid  string. 

In  determining  whether  a  substitution  should  take  place  in  the  sub¬ 
programs,  blanks  are  ignored.  That  is,  all  statements  in  the  input  pro¬ 
grams  are  compressed  as  though  no  blanks  were  left  anywhere.  Of  course, 
blanks  are  not  deleted  from  hollerith  fields. 
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2.3.  EXAMPLES 

1.  A=B(l,I)+SQRT(D(2,2))-F(I,J+l,K)+2  *G 

2.  IF(Z.LT. .3)  CALL  EXIT 

In  the  above  two  statements,  only  the  following  are  valid  strings: 


* 

A 

* 

F 

EvM) 

* 

I 

* 

B 

J+l 

* 

I 

* 

J 

SQRT(D(2,2)) 

* 

K 

* 

SORT 

* 

G 

* 

D(2,2) 

* 

Z 

* 

D 

* 

EXIT 

F(I,  J+l  ,K) 

Two  passes  are  normally  made  through  each  statement  in  a  pro¬ 
gram.  A  statement, including  all  its  continuation  cards,  is  read,  and  blanks 
are  deleted.  During  the  first  pass  every  valid  string  is  checked  for  sub¬ 
stitution,  and  the  substitution,  if  any,  is  made.  During  the  second  pass  all 
identifiers  and  all  variables  having  a  completely  numerical  subscript  are 
written  on  tape  along  with  the  new  source  line  in  which  each  variable 
appears.  Those  valid  strings  listed  above  with  asterisks  preceding  them 
would  have  been  recorded  on  tape.  Not  only  are  the  name  and  source  line 
of  each  identifier  and  numerically  subscripted  variable  recorded  on  tape, 
but  the  fact  that  it  appeared: 

1.  In  a  type  statement  or  specification  statement  as  (1)  complex, 
(2)  integer,  (3)  external,  (4)  logical,  (5)  real,  (6)  double  pre¬ 
cision,  (7)  equivalence,  (8)  parameter,  (9)  data. 

2.  On  the  left  side  of  an  equal  sign. 

3.  As  an  argument  of  (1)  the  subprogram  itself,  (2)  a  referenced 
subroutine,  (3)  a  referenced  function. 

4.  As  the  name  of  the  subprogram  in  a  subroutine  statement  or 
function  statement. 
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5.  As  the  name  of  a  referenced  subroutine  or  referenced  function. 

6.  In  a  statement  which  specifies  a  variable  to  be  dimensioned 
(in  a  common  or  type  statement  with  a  numerical  subscript  or 
in  a  dimension  statement),  in  which  case  its  dimension  and  the 
number  of  characters  in  its  dimension  are  also  recorded. 

7.  In  a  common  statement  as  either  (1)  a  block  common  name  or 
(2)  a  variable  name,  in  which  case  its  common  block  or  blank 
is  also  recorded. 

8.  As  a  variable  with  a  completely  numerical  subscript,  in  which 
case  its  subscript  and  the  number  of  characters  in  its  subscript 
are  also  recorded. 

9.  In  a  read  statement. 

10.  Simply  as  a  variable  in  none  of  the  above. 

This  information  is  recorded  at  the  time  of  appearance  for  each  and  every 
appearance  of  every  identifier  and  numerically  subscripted  variable. 

Pass  two  of  Program  One  edits  each  statement  after  substitutions,  if 
any  are  made,  in  the  following  manner;  Single  blanks  are  inserted  before 
and  after  FORTRAN  logical  operators  and  the  connectors  =,  +,  and  -,  fol¬ 
lowing  a  comma  and  neither  before  nor  after  the  connectors  *,  /, 

Blanks  do  not  occur  in  column  7  except  for  continuation  cards,  when  they 
always  occur.  Identifier  names  are  never  split  between  one  card  and  its 
continuation  card.  The  continuation  card  mark  in  column  6  is  always  an 
apostrophe.  The  following  statements  are  not  edited  except  that  the  con¬ 
tinuation  marks  are  changed  to  apostrophes  and  information  is  placed  in 
columns  73-80;  comment  cards,  FORMAT  statements,  CONTINUE  cards, 
and  END  statements.  The  compound  or  logical  IF  statement  is  split  so  that 
the  second  half  of  the  compound  statement  begins  on  a  continuation  card. 

On  all  cards,  columns  73-76  will  contain  the  first  four  characters  of  the 
subprogram  name  appearing  as  the  first  name  on  the  FOR  card  for  each 
subprogram.  Columns  77-80  will  contain  a  4  digit  sequence  number  start¬ 
ing  from  1  for  each  subprogram  and  increasing  by  1. 
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A  new  sequenced  source  deck  will  be  punched.  A  listing  of  the  source 
deck  will  be  printed  with  the  additional  information  that  the  first  line  of 
each  statement  where  a  substitution  has  been  made  will  be  so  noted.  Any 
statement  which  required  more  than  19  continuation  cards  after  editing  will 
be  so  noted,  and  the  additional  cards  will  be  punched  and  listed.  All  strings 
for  which  substitutions  were  requested  that  did  not  appear  in  the  subprogram 
just  piocessed  will  be  noted  at  the  end  of  the  listing  for  that  subprogram. 

In  addition  to  the  pertinent  description  of  Program  One  given  above, 
the  following  unrelated  facts  concerning  this  program  must  be  stated  and 
understood  if  it  is  to  be  used  properly: 

1.  Statement  numbers  are  in  no  way  affected  by  the  program. 

2.  This  program  does  not  have  extensive  checking  features,  and 
a  non-compilable  statement  may  create  an  endless  loop. 

3.  Each  subprogram  must  begin  with  a  FOR  card  with  column  one 
blank  (the  7-8  punch  removed)  and  end  with  an  END  statement. 

4.  Superfluous  real,  integer,  and  equivalence  statements,  created 
by  virtue  of  variable  substitutions,  will  not  be  eliminated. 

5.  No  new  real-  or  integer-type  statements  will  be  generated  to 
correct  for  an  unintentional  change  of  type  when  a  name  sub¬ 
stitution  changes  the  type. 

6.  A  name  substitution  will  take  place  whenever  that  name  appeal  s 
regardless  of  whether  it  is  subscripted. 

7.  A  variable  plus  subscript  substitution  will  take  place  only  when 
the  variable  and  subscript  agree  exactly.  However,  no  sub¬ 
stitution  will  take  place  if  the  subscript  is  numeric  and  appears 
in  a  dimension  statement  or  any  statement  which  specifies  its 
dimensionality.  Thus,  the  replacement  of  R(152)  by  PCYCLE 
would  not  occur  in  the  staterr-nt  COMMON  R(152). 

8.  The  appearance  of  any  variable  followed  by  a  left  parenthesis 
not  appearing  in  a  specification  statement  defining  its  dimen¬ 
sionality  will  be  checked  to  determine  whether  it  has  been 
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previously  dimensioned.  If  not,  it  will  be  assumed  to  be  a 
function  reference.  Thus,  on  substituting  J(5)  for  B,  if  C  has 
not  been  previously  dimensioned,  it  will  be  considered  a  function 
and  will  not  be  otherwise  rote  i  as  a  possible  source  error. 

9.  This  program  cannot  accept  comment  cards  embedded  between 
the  first  card  of  a  statement  and  the  last  continuation  card  of 
the  statement  Such  cards  must  be  removed  from  a  deck  before 
this  program  is  used.  Comment  cards  may  appear  any  place 
else  between  a  FOR  card  and  an  END  card  of  a  subprogram. 

10.  Anything  appearing  in  columns  73-80  on  comment  cards  is 
deleted,  and  the  sequencing  scheme  previously  described  is 
substituted. 

11.  The  number  of  +  signs,  -  signs,  and  logical  operators  appear¬ 
ing  in  a  compound  or  logical  IF  statement  must  not  exceed  20 
if  this  statement  is  to  bt  reproduced  correctly. 

Program  Two  processes  the  output  of  Program  One.  It  will  print 
out  the  following  information  for  each  subprogram: 

1.  The  subprogram  name  plus  an  18  character  tape  identification. 
This  identification  is  described  further  in  Section  2.4  of  this 
report. 

2.  An  alphabetic  listing  of  all  common  blocks  starting  with  blank 
common.  Fur  each  common  block  the  following  information  is 
printed. 

a.  A  listing  of  the  variables  within  each  block  in  order  of 
their  appearance  in  the  block,  and  for  each  variable  its 
octal  and  decimal  origin  relative  to  zero  within  the  block. 

b.  An  alphabetic  listing  of  the  various  within  each  block, 
and  for  each  variable  its  octal  and  decimal  origin  rela¬ 
tive  to  zero  within  the  block. 
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3.  An  alphabetic  listing  of  the  subroutines  and  fractions  referenced 

* 

within  this  subprogram  and  the  source  lines  in  which  they  were 
referenced. 

4.  An  alphabetic  listing  of  each  variable  appearing  in  the  sub¬ 
program.  For  each  variable  the  following  information  is  printed: 

a.  A  single  line  containing  the  variable  name,  the  block 
name  in  which  it  may  appear,  its  octal  and  decimal  origin 
relative  to  zero  within  that  block,  its  dimensions  (if  any), 
the  size  of  its  dimensions,  and  a  series  of  code  letters 
describing  its  "uses"  within  the  program  (i.e.  ,  its  type, 
whether  it  arrears  anywhere  within  the  subprogram  on 
the  left  side  of  an  equal  sign,  whether  it  appears  any¬ 
where  within  the  subprogram  as  the  argument  of  a  call, 
whether  it  appears  anywhere  in  the  subprogram  in  a  read 
statement,  etc. )  A  full  description  and  listing  of  these 
code  letters  appears  in  Section  2.4  of  this  report  and  on 
the  first  page  of  the  printed  output  of  Program  Two. 

b.  An  ordered  listing  of  every  source  line  in  which  the 
variable  appears  is  printed.  Following  each  source  line 
there  may  appear  the  same  coded  letters  as  are  des¬ 
cribed  in  (4a),  above,  if  the  variable  can  be  modified 
within  that  line  (by  appearing  on  the  left  side  of  an  equal 
sign  or  as  an  argument  of  a  function,  in  a  read  state¬ 
ment,  etc.). 

c.  For  each  variable  that  is  numerically  subscripted  within 
the  program,  a  sort  order  listing  for  each  numeral  sub¬ 
script,  followed  by  a  list  of  all  source  lines  in  which  it 
appears  with  this  subscript. 

After  all  subprograms  have  been  processed  and  each  of  the  above 
listings  has  been  printed,  the  program  prints  an  alphabetic  listing  of  all 
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3.  An  alphabetic  listing  of  the  subroutines  and  functions  referenced 
within  this  subprogram  and  the  source  lines  in  which  they  were 
referenced. 

4.  An  alphabetic  listing  of  each  variable  appearing  in  the  sub¬ 
program.  For  each  variable  the  following  information  j.c  printed: 

a.  A  single  line  containing  the  va  iuble  name,  the  block 
name  in  which  it  may  appear,  its  octal  and  decimal  origin 
relative  to  zero  within  that  block,  its  dimensions  (if  any), 
the  size  of  its  dimensions,  and  a  series  of  code  letters 
describing  its  ''uses"  within  the  program  (i.  e.  ,  its  type, 
whether  it  appears  anywhere  within  the  subprogram  on 
ike  left  side  of  an  equal  sign,  whether  it  appears  any¬ 
where  within  the  subprogram  as  the  argument  of  a  call, 
whether  it  appears  anywhere  in  the  subprogram  in  a  read 
statement,  etc. )  A  full  description  and  listing  of  these 
code  letters  appears  in  Section  2.4  of  this  report  and  on 
the  first  page  of  the  printed  output  of  Program  Two. 

b.  An  ordered  listing  of  every  source  line  in  which  the 
variable  appears  is  printed.  Following  each  source  line 
there  may  appear  the  same  coded  letters  as  are  des¬ 
cribed  in  (4a),  above,  if  the  variable  can  be  modified 
within  that  line  (by  apDearing  on  the  left  side  of  an  equal 
sign  or  as  an  argument  of  a  function,  in  a  read  state¬ 
ment,  etc. )  . 

c.  For  each  variable  that  is  numerically  subscripted  within 
the  program,  a  s«>rt  order  listing  for  each  numeral  sub¬ 
script,  followed  by  a  list  of  all  source  lines  in  which  it 
appears  with  this  subscript. 

After  all  subprograms  have  been  processed  and  each  of  the  above 
listings  has  been  printed,  the  program  prints  an  alphabetic  lifting  of  ail 
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identifiers  appe  mg  in  any  of  the  subprograms,  and  for  each  identifier, 
an  alphabetic  listing  of  the  subprograms  in  wh.  1  it  appeared.  The  sub¬ 
programs  are  listed  by  subroutine  name  (not  by  the  name  on  the  FOR 
card).  This  ends  the  output  description  of  Program  Two. 

Both  Program  One  and  Program  Two  use  logical  unit  10  to  store  and 
process  the  identifier  information.  A  complete  description  of  the  format 
and  information  stored  on  logical  unit  10  by  Program  One  is  contained  in 
Section  2.  5  of  this  report.  Table  lengths  for  Program  Two  are  also  given 

in  that  section. 

2.4.  INPUT 

It  is  suggested  that  FORTRAN  code  sheets  be  used  for  all  data  input, 
since  the  data  format  has  been  so  designed.  There  are  three  types  of  data 
input  for  Program  One.  Data  type  1  is  a  single  card  having  two  fields. 

It  precedes  a  complete  run  of  Program  One  which  may  consist  of  one  or 
more  subprograms.  Field  1  -  Column  1  (format  II)  controls  logical  unit  10; 

if  field  1  is  blank,  logical  unit  10  will  be  rewound.  This  implies  either 
that  the  unit  is  blank  or  that  information  previously  stored  on  this  unit  is 
no  longer  desired.  If  field  1  contains  a  one,  it  is  assumed  that  logical 
unit  10  had  information  previously  stored  on  it  from  a  previous  run  of 
Program  One;  then  logical  unit  10  will  be  positioned  forward  to  the  end  of 
the  previously  stored  information,  ready  to  receive  additional  information. 
Field  1  -  Columns  2-6  inclusive  are  blank.  Field  2  -  Columns  7-24  (for¬ 
mat  3A6)  is  the  secondary  identification  valid  for  all  subprograms  of  this 
run.  It  may  contain  any  information  desired,  e.g. ,  the  date.  It  is  in¬ 
tended  that  a  program  will  be  written  to  edit  logical  unit  10,  deleting  one 
version  of  a  subprogram  and  accepting  another  version  of  this  same  sub¬ 
program  by  utilizing  the  uniqueness  of  the  18  character  identification  in 
field  2  between  runs  of  Program  One. 
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For  Data  type  2,  one  or  more  cards  must  precede  each  subprogram, 
but  no  subprogram  will  require  more  than  15  of  these  cards.  Fifteen 
"blocks"  have  been  set  aside,  internal  to  Program  One,  for  "substitutions." 
(A  "substitution"  is  defined  as  the  combination  of  a  source  program  string 
to  be  changed  and  its  replacement  string. )  A  block  in  which  no  substitu¬ 
tion  has  ever  been  placed  or  one  in  which  a  substitution  has  been  deleted 
without  replacement  is  considered  void.  A  block  that  is  not  void  is  full. 

To  "delete"  means  to  wipe  out  a  substitution.  It  can  only  be  brought  back 
by  fully  restating  the  substitution  on  a  new  type  2  card,  as  described 
below.  A  type  2  card  may  serve  any  of  the  following  purposes: 

1.  It  may  fill  a  void  block  with  a  substitution.  The  block  is 
then  full.  The  substitution  then  applies  through  all  sub¬ 
programs  that  follow  it  unless  deleted  or  replaced  by  an¬ 
other  type  2  card. 

2.  It  m-y  delete  a  particular  full  block.  Tnat  block  is  then  void. 

3.  It  may  delete  all  existing  full  blocks.  All  blocks  will  then  be 
void. 

4.  It  may  delete  a  substitution  from  a  full  block  and  replace  it 
with  another  substitution  in  the  same  block.  That  block  is 
then  still  full.  The  new  substitution  will  then  apply  through 
all  subprograms  that  follow  it  unless  it  is  specifically  deleted 
or  replaced  by  still  another  type  2  card. 

Type  2  cards  consist  of  three  fields,  hereafter  referred  to  as  FI, 

F2,  and  F3.  FI  (Format  A6)  delineates  the  last  data  type  2  card  by 
having  some  uon-blank  character  in  any  of  columns  1-6  of  the  data  card, 
(except  that  column  1  must  not  have  a  7-8  punch).  That  is,  columns  1-6 
must  be  blank  for  all  but  the  last  data  type  2  card  preceding  each  sub¬ 
program.  4 
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F2  (format  71A1)  is  either  entirely  blank,  in  which  case  it  sets  the 
block  specified  in  F3  void,  or  contains  the  substitution  in  columns  7-77, 
inclusive.  If  non-blank.,  the  general  form  of  F2  is  SOURCE  PROGRAM 
STRING  TO  BE  MODIFIED  =  ITS  REPLACEMENT  STRING.  The  following 
rules  must  be  observed  in  the  specification  of  F2. 

1.  F2  must  start  in  column  7.  (It  is  considered  blank  if  column 
7  of  F2  is  blank. ) 

2.  The  two  strings  must  be  separated  by  an  equal  sign. 

3.  There  must  not  be  any  blanks  between  column  7  and  the  last 
character  of  the  replacement  string. 

4.  Column  78  must  be  blank. 

For  F3  (Format  12),  columns  79  and  80  must  contain  blanks  or  the 
numbers  01-15  or  the  number  99.  F3  may  be  blank  if  and  only  if  FI  is 
non-blank  and  F2  is  blank,  i.e.  ,  to  express  the  end  of  type  2  data  without 
supplying  any  information  on  F2  and  F3.  In  general,  F3  contains  the  block 
numbers  01-15  for  which  an  F2  substitution  or  deletion  or  replacement 
substitution  is  being  made.  The  information  for  a  block  need  not  be 
ordered,  i.e. ,  information  for  block  07  may  be  supplied  before  information 
for  block  01,  and  blocks  may  be  arbitrarily  full  or  void.  If  F3  contains  a  99, 

all  15  blocks  are  void.  F2  would  be  blank  if  F3  contained  a  99,  and  FI 
might  or  might  not  be  blank.  It  is  perfectly  valid  for  additional  type  2 
cards  to  follow  an  F3  containing  a  99  (in  which  case  FI  would  be  blank 
when  F3  contained  a  99). 

If  all  15  blocks  are  void  upon  return  from  the  program  that  processes 
the  type  2  cards,  then  Program  One  will  bypass  its  first  pass  (the  substi¬ 
tution  phase)  and  will  presumably  run  twice  as  fast.  This  implies  that  all 
one  is  interested  in  is  the  cross  reference  dictionary  generated  in  pass 
two.  If  this  is  so,  then  the  data  type  2  cards  preceding  each  subprogram 
should  consist  of  a  single  card  with  perhaps  a  single  character  in  column  1. 
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(However,  the  single  character  must  not  be  a  7-8  punch,  which  transfers 
control  to  the  1108  operating  system.) 

Data  type  3  cards  are  subprogram  source  cards  to  be  processed. 

They  must  start  with  a  FOR  card  with  the  column  1  blank  (that  is  the  7-8 
punch  must  be  removed).  Following  column  1  there  must  be  at  least  one 
blank  preceding  and  at  least  one  blank  following  the  "FOR"  and  before  its 
Deck  Name.  Just  prior  to  the  FOR  card  there  must  be  a  data  type  2  card 
with  at  least  one  character  other  than  a  blank  in  columns  1-6  (and  no';  a 
7-8  punch  in  column  1).  Following  the  modified  FOR  card  are  the  so  iree 
cards,  which  must  be  compilable.  The  source  deck  must  end  with  an  END 
card. 

Program  Two  requires  no  card  input.  A  request  for  logical  unit  10 
(ASG)  must  be  supplied  by  the  user  for  both  Program  One  and  Program 
Two. 

2.  5.  APPENDIX:  PROGRAM  DETAILS  ■ 

Table  I  contains  information  describing  the  code  letters  which  define 
the  "use"  of  an  identifier  (see  steps  (4a)  and  (4b)  of  the  Program  Two  out¬ 
put  description  in  Section  2.2),  and  the  code  numbers  written  on  tape  by 
Program  One  which  identify  that  "use.  " 

Any  combination  of  the  code  letters  appearing  in  column  1  of  Table  I 
may  be  printed  in  the  Program  Two  output  previously  described  in  step  (4a) 
of  Section  2.3,  but  only  those  code  letters  in  Table  I  which  have  an  asterisk 
preceding  them  will  be  printed  out  as  described  in  step  (4b).  The  latter 
code  letters  reference  only  those  statements  which  can  modify  a  variable. 

The  information  written  by  Program  One  on  logical  unit  10  is  written 
in  unformatted  I/O  (equivalent  to  binary)  and  consists  of  six  words  per 
record.  The  six  words  make  up  four  fields.  Field  1  is  in  field  data  and  is 
either  (1)  the  program  name  appearing  on  the  FOR  card,  as  the  first  word 
of  the  first  and  last  record  of  each  subprogram,  or  (2)  an  identifier  appear- 
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Table  I 

DESCRIPTION  OF  CODE  LETTERS 


Code 

Letter 

Code 

Number(k) 

U 

-  89 

*  D 

-  90 

-  9  i 

7. 

-  92 

*  X 

-  93 

*  H 

-  94 

-  93 

*  S 

-  96 

B 

-  97 

*  R 

-  98 

-  99 

E 

-100 

C 

-iOl 

i 

-102 

■*  o 

-103 

L 

-104 

F 

-  i05 

P 

-106 

A 

-107 

*x  G 

-ioe 

-i8*k<0 


0  <k  5i8 


-2i8^k<-200 


0 


De«c  ripl  ion 

Variable  Appears  in  a  parameter  statement 
Variable  appears  in  a  data  statement 
Name  appears  on  FOR  Card 

Variable  appears  in  a  complex  type  statement 
Variable  appears  on  left  side  of  an  equal  sign 
Variable  appears  in  a  common  block 
Block  name  appears  in  a  common  statement 
Variable  appears  as  an  argument  of  a  cali 
Variable  appears  in  a  double  precision  type  statement 
Variable  appears  in  a  read  statement 

Subroutine  name  on  subroutine  or  function  statement  card 
Variable  appears  lr  an  external  type  statement 
Subroutine  nar.ie  appears  in  a  call  statement 
Variable  name  appears  in  an  integer  type  statement 
Variable  name  appears  in  an  equivalence  type  statement 
Variable  appeals  in  a  logical  type  statement 

Function  name  appear  in  an  arithmetic  statement  or  IF  statement,  or 
on  either  sids  of  the  equal  sign  in  an  arithmetic  function  statement 
(but  not  on  a  function  statement  card  described  in  -99  above) 

Variable  appears  in  a  real  type  statement 

Variable  appears  as  a  program  argument  on  a  function  or  subroutine 
statement  card 

Variable  appears  as  the  argument  of  a  referenced  function  as  des¬ 
cribed  in  -105  above 

Variable  appears  in  a  dimensior.  statement  or  in  a  statement  describing 
its  dimensionality  where  |k|  is  the  number  of  characters  in  its  dimen¬ 
sion  description  including  parentheses,  commas,  and  digits 

Variable  appears  in  a  statement,  other  than  one  describing  its  dimen¬ 
sionality  or  a.'>  equivalence  statement,  with  a  numerical  subscript 
where  k  is  the  number  of  characters  in  its  subscript  including  both 
parentheses,  commas,  and  digits 

Variable  appears  in  an  equivalence  statement  with  a  subscript  where 
I k-200 1  is  the  number  of  charac'.ers  in  its  subscript,  including  both 
parentheses,  commas,  and  digits 

Not  specified  above  (as  the  variables  A  or  B  in  the  arithmetic  state¬ 
ment  Y  =A  +  B) 
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ing  in  the  program.  Field  1  is  left  adjusted  and  filled  with  blanks  on  the 
right  if  the  identifier  is  less  than  six  characters. 

Field  2  is  in  field  data  and  consists  of  three  words,  either: 

1.  The  18  character  identification  specified  on  the  data  type  1 
card  of  a  specific  run  of  Program  One.  It  will  appear  in  the 
second  field  of  the  first  record  of  every  subprogram  processed 
during  this  run. 

2.  A  variable  dimension,  left  adjusted,  with  no  embedded  blanks, 
containing  both  parentheses  and  their  contents,  and  filled  with 
blanks  on  the  right.  Field  2  will  contain  this  information  only 
if  the  variable  in  field  1  is  appearing  in  a  statement  specifying 
its  dimensionality  at  the  time  this  record  is  written. 

3.  A  variable  numerical  subscript,  left  adjusted,  with  no  embedded 
blanks  (as  described  in  step  2  above).  Field  2  will  contain  this 
information  only  if  the  variable  in  field  1  is  appearing  in  any 
statement  other  than  one  specifying  its  dimensionality  at  the 
time  this  record  is  written  and  is  numerically  subscripted. 

4.  A  block  name  (which  may  be  blank  for  blank  common).  Word  1 
of  field  2  will  contain  this  information  only  if  the  variable  is 
appearing  in  a  common  statement  at  the  time  the  record  is 
written;  then  Words  2  and  3  of  field  2  will  be  blank.  If  field  1 
contains  the  block  name,  the  first  word  of  field  2  will  also 
contain  that  block  name.  All  three  words  of  field  2  contain 
blanks. 

Field  3,  the  fifth  word  of  the  record,  is  the  code  number  appearing 
in  the  second  column  of  Table  I.  The  last  record  of  each  subprogram  con¬ 
tains  a  -200  in  field  3  (generated  by  the  END  statement).  Field  4  contains 
the  source  line  number  in  which  the  identifier  in  field  1  appears.  Following 
the  last  record  of  the  last  subprogram  is  one  additional  six  word  record 
containing  a  -200  in  its  fifth  word.  There  is  no  end-of-file  following  this. 
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One  can,  therefore,  determine  the  end  of  information  by  seeking  two  con¬ 
secutive  records  each  of  which  contains  -200  in  its  fifth  word. 

If  a  variable  appears  in  a  statement  specifying  its  dimensionality  or 
is  numerically  subscripted,  it  will  appear  in  two  consecutive  records. 

For  example,  the  numerically  subscripted  variable  A(10)  in  the  statement, 
A(10)=Y  appearing  in  line  17  of  a  subprogram  appear  s  in  two  records  as: 


Field  1 

Field  2 

Field  3 

Field  4 

Record  n 

A 

U0) 

4 

17 

Record  n+1 

A 

-93 

17 

Record  n+2 

Y 

0 

17 

The  purpose  of  the  above  description  is  to  provide  the  user  with  an 
opportunity  to  write  his  own  specialized  cross  reference  dictionary  (i.e.  , 
his  own  variation  of  Program  Two.  For  example,  he  may  wish  to  print  out 
all  variables  appearing  in  READ  statements  throughout  the  program  in  the 
order  in  which  they  are  read. 

Certain  assumptions  had  to  be  made  in  Program  Two  concerning  the 
internal  allocation  of  space  for  the  many  tables  required  by  the  cross  ref¬ 
erence  dictionary.  This  allocation  i"  further  complicated  by  a  possible 
variation  in  core  space  that  will  be  available  to  the  user  in  the  near  future. 

It  is  currently  assumed  that  the  total  data  space  must  not  exceed  40,  000. 
That  is,  the  sum  of  the  "total  space"  created  by  the  parameter  occurrences 
listed  in  Table  II  and  fixed  working  memory  must  not  exceed  40,  000, 
However,  as  long  as  the  total  does  not  exceed  40,  000,  any  one  or  more  of 
the  five  parameters  may  be  modified  within  the  program  in  the  event  the 
individual  tables  to  which  they  refer  are  exceeded.  This  requires  the 
change  of  a  single  PARAMETER  card  in  Program  Two,  which  states: 
PARAMETER  Zl  =  15000,  Z2  =  500,  Z3  =  1000,  Z4=2500,  Z5=1000,  where 
Z1-Z5  and  the  number  of  appearances  of  Z1-Z5  are  as  defined  in  Table  II. 
The  "limit"  column  in  Table  II  refers  to  the  maximum  appearances  or  refer¬ 
ences  or  number  of  such  items  in  any  one  subprogram. 
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Para¬ 

meter 

Z1 

Z2 

Z3 

Z4 

Z  5 


Table  H 

DEFINITION  OF  ENTRIES  IN  PARAMETER  CARD 


Description 

Number  of  references  of  vari¬ 
ables  with  or  without  sub¬ 
scripts 

Number  of  unique  common 
block  names 

Number  of  appearances  of 
variables  with  unique  sub- 
scfipts 

Number  of  appearances  of 
unique  identifiers 

Number  of  appearances  of 
unique  variables  that  are 
dimensioned 


Number 
of  tables 
requiring 


Limit 

this  para¬ 
meter  = 

Total 

space 

15000 

x  1 

15, 000 

500 

x  1 

500 

1000 

x  6  = 

6,  000 

2500 

x 

rf* 

ii 

10,000 

1000 

x  5 

5,000 

36,  500 


Fixed  working  memory  (not 

subject  to  change)  3’  500 

40,000 
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SECTION  III 

DATA-LINK  EDITING  ROUTINES 

3.  1.  INTRODUCTION 

In  order  to  facilitate  use  of  the  computer  codes  developed  in  the 
NEIS  on  the  CDC-6600  computer  at  the  Air  Force  Weapons  Laboratory,  a 
data  link  was  installed  between  the  SNEL  facility  in  San  Diego  and  the  AFWL 
facilities  in  Albuquerque.  The  link  terminals  include  an  IBM-1978  reader- 
printer  at  SNEL,  which  communicates  with  an  IBM-7702  tape  terminal  at 
AFWL.  The  format  of  the  information  which  can  be  transmitted  or  received 
by  the  1978  is  unfortunately  incompatible  with  the  requirements  of  the 
Chippewa  operating  system  for  the  6600  computer,  so  editing  procedures 
are  required  as  an  intermediate  step  in  the  data  transmission  process. 

Two  programs,  called  TOGA  and  FROGA,  were  written  for  the  AFWL 
CDC-16QA  computer  to  perform  these  editing  functions. 

3.2.  TOGA 

The  TOGA  program  edits  BCD  output  tapes  prepared  by  the  6600, 
producing  modified  tapes  which  can  be  read  by  the  7702,  and  printed  on  the 
1978.  The  following  editing  functions  are  accomplished: 

1.  Printer  carriage  control  characters  written  in  FORTRAN  form 
are  blank,  0,  and  1,  indicating  single  space,  double  space, 
triple  space,  and  page  eject.  These  actions  precede  the  printing 
of  the  line  in  which  they  occur.  The  corresponding  1978  control 
characters  are  /  ,  S,  T  and  A,  respectively,  with  the  difference 
that  the  actions  follow  the  printing  of  the  line  in  which  they  occur. 
The  TOGA  code  therefore  r  eplaces  each  FORTRAN  carriage 
control  character  by  its  1978  counterpart,  which  is  then  placed 
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at  the  head  of  the  preceding  line.  If  the  first  character  in  the 
line  is  already  a  1978  character,  no  action  is  taken;  if  it  is  an 
illegal  character,  a  single  space  (/)  is  supplied. 

2.  The  1978  requires  each  line  to  be  terminated  by  a  segment 
mark,  which  is  a  special  BCD  character.  This  character  is 
placed  at  the  end  of  each  line,  after  any  trailing  blanks  are 
deleted. 

3.  The  physical  records  transmitted  over  the  link  are  blocked  to 
some  extent.  The  number  of  lines  per  block  varies  from  two 
to  seven,  depending  upon  the  total  number  of  characters,  which 
may  not  exceed  328. 

4.  The  length  of  a  line  printable  on  the  1978  is  120  characters, 
plus  the  leading  control  character  and  the  trailing  segment 
mark.  If  the  input  lines  exceed  this  length,  the  overflow  goes 
onto  a  second  line  flagged  by  a  leading  exclamation  point,  follow¬ 
ing  which  the  carriage  is  double -spaced. 

5.  Input  records  less  than  three  characters  long  are  considered 
noise  and  are  skipped. 

3.3.  FROGA 

The  FROGA  program  edits  information  transmitted  by  the  1978  into  a 
form  acceptable  for  job  input  by  the  Chippewa  system.  The  following  func¬ 
tions  are  performed: 

1.  The  transmitted  records,  except  possibly  the  last,  consist  of 
four  80  character  card  images,  each  followed  by  a  segment 
mark.  These  records  are  read  one  at  a  time.  For  each  card 
image,  the  segment  mark  and  any  trailing  blanks  are  deleted, 
the  iatter  two  at  a  time. 

2.  The  card  images  are  moved  to  an  output  area.  The  end  of  each 
card  image  is  indicated  by  the  following  rule:  The  number  of 
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characters  for  each  card  image  must  be  a  multiple  of  10,  and 
the  last  two  of  these  must  be  6-bit  (binary)  zeros. 

As  the  BCD  characters  are  moved  to  the  output  area,  they  are 
translated  into  the  Chippewa  display  code  by  table  lookup. 

The  size  of  the  output  buffer  is  5120  characters.  The  buffer  is 
written  onto  an  output  "display  tape"  when  full,  or  when  an  end- 
of-record  flag  is  received.  The  end-of-record  flag  is  an  apos¬ 
trophe  in  column  1  of  a  card  image.  If  the  buffer  is  just  filled 
before  the  end-of-record  flag  is  received  (an  unlikely  occur¬ 
rence),  a  special  short  record  consisting  of  eight  binary  zeros 
(48  bits)  is  written.  Either  this  special  short  record  or  the 
normal  short  record  (less  than  512  ten  character  words)  in¬ 
dicates  to  the  Chippewa  system  that  a  "logical  record"  is 
complete. 

End-of-file  records  are  written  on  the  display  tape  when  a 
second  apostrophe  is  received  following  the  first,  e.g.  ,  :u 
column  2  of  the  same  card  image.  This  end-of-file  separates 
one  job  from  another.  A  third  apostrophe  indicates  end  of 
transmission,  as  well  as  end-of^record  and  end-of-file. 

When  only  one  job  is  transmitted  at  a  time,  only  single  oi 
triple  apostrophe  flags  are  used  (the  preferrsd  procedure). 

On  end-of-transmis  sion,  two  end-of-file  records  are  written 
on  the  display  tape. 
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column  2  of  the  same  card  image.  This  end-of-file  separates 
one  job  from  another.  A  third  apostrophe  indicates  end  of 
transmission,  as  well  as  end-of-record  and  end-of-file. 

When  only  one  job  is  transmitted  at  a  time,  only  single  or 
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